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COSMIC  RAY  EFFECTS  ON  MICROELECTRONICS, 

PART  I:  THE  NEAR-EARTH  PARTICLE  ENVIRONMENT 

1.0  Introduction 

A  number  of  ground-based  experiments  have  recently  shown  that  a  single 
intensely-ionizing  particle  can  change  the  logic  state  of  modern  electronic 
circuits  of  the  kind  used  as  memories  on  satellites.  Soft  errors  (also 
called  soft  upsets  or  single  event  upsets)  have  been  observed  on  more  than  a 
dozen  satellites  .  The  soft  errors  on  three  of  these  satellites  have  been 
conclusively  attributed  to  single  intensely-ionizing  particles.  Besides 
soft  errors,  single  intensely  Ionizing  particles  have  been  shown,  in 
laboratory  experiments,  to  cause  latchup  and  even  to  do  permanent  damage  to 
the  microcircuits. 

Intensely-ionizing  particles  may  be  produced  locally,  in  the  electronic 
device  Itself,  as  a  product  of  nuclear  reactions  or  they  may  come  directly 
from  outside  the  spacecraft.  Ever  when  the  Intensely-Ionizing  particle  is 
the  product  of  a  nuclear  reaction,  that  reaction  is  usually  initiated  by  a 
more  lightly  ionizing  particle  that  came  from  outside  the  spacecraft. 

The  objective  of  this  study  is  to  begin  addressing  this  problem  by 
developing  the  tools  needed  to  estimate  the  rate  at  which  soft  errors  can  be 
expected  to  occur  on  various  spacecraft  exposed  to  the  natural  space 
environment.  The  first  step  is  to  develop  a  model  of  the  energetic  particle 
environment  near  earth  that  is  accurate  and  yet  easy  to  use. 

This  report  will  be  followed  by  additional  reports.  One  will  describe 
the  way  in  which  the  earth's  magnetic  field  has  modulated  the  energy  spectra 
of  particles  reaching  any  satellite.  A  second  report  will  describe  how  the 
energy  spectra  and  elemental  composition  of  these  particles  are  modified  in 
passing  through  spacecraft  walls  to  reach  the  electronic  components  inside. 
The  results  of  these  three  reports  can  then  be  combined  with  measured  or 
estimated  operational  cross  sections  for  the  various  single-particle  effects 
on  microelectronics  to  compute  their  expected  rates  on  various  spacecraft. 

This  report  describes  simple  analytic  models  for  the  energy  spectra  and 
elemental  compositions  of  the  various  components  of  ionizing  particle 
radiation  in  the  vicinity  of  the  earth  that  are  as  accurate  as  the  data  will 
allow.  The  models  are  based  on  an  exhaustive  review  of  the  available  data. 
From  the  length  of  this  report,  it  can  be  seen  that  a  substantial  data  base 
exists  on  the  energetic  particle  environment.  Even  so  there  are 
deficiencies  In  the  data  base  required  to  accurately  estimate  the  rates  of 
single  particle  effects. 

This  situation  has  led  us  to  adopt  the  following  philosophy  in  modeling 
the  environment.  Where  the  data  base  is  adequate,  the  model  gives 
"most-probable"  spectra  and  compositions.  When  the  component  is  variable,  a 
worst  case,  at  a  90  per  cent  confidence  level,  is  given.  In  those  cases 
where  the  data  base  is  inadequate ,  we  can  only  speculate  what  the  conditions 
might  be.  Such  speculation  would  lead  us  to  construct  credible  worst-case 
models  that  are  quite  severe  and  therefore  pessimistic  from  the  spacecraft 
designers  point  of  view.  To  avoid  provoking  undue  expense  in  spacecraft 
design,  we  have  adopted  an  optimistic  philosophy.  In  cases  where  the  data 
base  is  Inadequate,  we  have  modeled  whatever  data  actually  exist,  ignoring 
the  untested  possibilities.  This  guarantees  the  user  that  his 
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spacecraft  will  actually  experience  an  environment  as  severe  as  the  one 
described  here.  The  design  effort  expended  by  using  this  model  will  then 
not  have  been  wasted.  It  Is,  of  course,  possible  that  some  of  the  untested 
speculations  may  prove  correct,  leading  to  a  far  more  hostile  environment 
than  described  here.  Until  the  necessary  experiments  can  be  carried  out, 
the  spacecraft  designer  must  simply  take  some  risks. 

Spacecraft  operating  near  earth  may  be  bombarded  by  energetic  charged 
particles  that  are  trapped  In  the  earth's  radiation  belts.  Spacecraft  may 
also  be  bombarded  by  cosmic  rays,  particles  from  solar  flares  and  particles 
accelerated  In  the  Interplanetary  medium,  all  of  which  come  from  great 
distances. 

The  contribution  each  component  makes  to  the  total  particle  flux 
bombarding  a  spacecraft  is  complicated  by  the  presence  of  the  earth's 
magnetosphere  shown  In  Figure  1.1.  The  Intensity,  energy  and  elemental 
composition  of  the  trapped  radiation  varies  enormously  with  position  In  the 
radiation  belts.  To  reach  a  spacecraft  Inside  the  magnetosphere,  particles 
coming  from  great  distances  must  penetrate  the  earth's  magnetic  field. 

Their  ability  to  do  so  depends  on  their  momentum  divided  by  their  electrical 
charge.  The  larger  this  ratio,  the  deeper  they  can  penetrate. 

In  the  mooels  presented  here,  we  will  describe  the  trapped  radiation  as 
It  Is  found  In  the  radiation  belts.  The  cosmic  rays,  solar  flare  particles 
and  particles  from  the  Interplanetary  medium  will  be  described  as  they  are 
found  outside  the  magnetosphere  In  the  Interplanetary  medium  near  the  orbit 
of  the  earth.  A  later  report  will  describe  how  these  components  are 
modulated  prior  to  reaching  the  orbit  of  a  satellite  In  the  magnetosphere. 

As  pointed  out  In  the  beginning  paragraphs  of  this  Introduction,  It  Is 
the  Intensely-Ionizing  particles  that  cause  single  particle  effects  on 
microelectronics.  The  Intensity  with  which  a  charged  energetic  particle 
ionizes  matter  varies  approximately  as  the  square  of  the  particle's 
electrical  charge  divided  by  the  square  of  Its  velocity.  When  a  particle  Is 
Ionizing  Intensely  enough  to  produce  a  single  particle  effect  directly,  It 
will  be  far  more  effective  In  doing  so  than  a  particle  that  must  produce  a 
nearby  nuclear  reaction  with  an  Intensely-Ionizing  product.  This  difference 
in  effectiveness  Is  about  1C6 ,  so  the  energy  spectra  and  elemental 
compositions  of  energetic  particles  In  the  natural  environment  are  very 
Important  for  the  estimation  of  these  effects. 

The  energy  spectra  presented  here  are  differential  energy  spectra.  They 
give  the  particle  flux  per  unit  energy  as  a  function  cf  the  particle's 
energy.  The  units  of  energy  are  millions  of  electron  volts  per  atomic  mass 
unit  (MeV/u)  or  billions  of  electron  volts  per  atomic  mass  unit  (GeV/u)  (see 
Rossi,  1964,  Appendix  E  for  an  explanation  of  electron  volt).  This  way  of 
expressing  energy  is  useful  because  It  means  that  particles  with  the  same 
energy  also  have  the  same  velocity  regardless  of  their  atomic  mass.  Many  of 
the  properties  of  the  various  elemental  spectra  are  identical  when  this 
energy  scale  Is  used.  The  units  of  flux  are  particles  per  square  meter  - 
second  -  steradlan  -  MeV/u  (m2.sec.ster.MeV/u).  The  steradian  Is  a  unit 
of  solid  angle. 
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The  elemental  composition  of  most  of  the  energetic  particle  components 
Is  similar  to  the  universal  composition  of  matter  as  oetermined  from  the 
study  of  meteorites,  the  sun  and  the  stars.  Figure  1.2  shows  the  relative 
abunoances  of  the  elements  In  nature  (Cameron,  1980).  As  can  be  seen  In 
Figure  1.2  the  elements  are  -  93.6  per  cent,  hydrogen  ~  0.3  per  cent 
helium  and  -  0.14  per  cent  elements  carbon  and  heavier.  Iron  Is  about  E 
per  cent  of  oxygen  and  the  elements  beyond  nickel  are  very  rare.  This  is 
approximately  the  composition  seen  In  solar  flare  particles,  though  the 
actual  composition  varies  a  lot  from  flare  to  flare.  The  galactic  cosmic 
ray  composition  Is  qualitatively  similar  to  Figure  2.1,  but  differs 
considerably  in  detail.  The  compositions  of  particles  accelerated  In  the 
Interplanetary  meoium  and  trapped  In  the  magnetosphere  are  profoundly 
altered  by  special  physical  effects. 

For  those  who  do  not  have  a  background  In  space  science,  but  wish  to 
know  more  about  the  subject,  we  recommend,  “Cosmic  Rays"  by  Bruno  Rossi 
(1964),  "Space  Physics"  by  Steve  White  (1970),  and  "Introduction  to  Space 
Science"  by  Wllmont  Hess  and  Gilbert  Mead  ( 1960 ) . 

For  users  of  this  report  who  are  Interested  only  In  the  model  Itself, 
the  details  have  been  collected  in  Appendix  1.  This  appendix  gives  all  the 
equations  required  to  compute  the  flux  levels  expected  under  various 
conditions  In  the  near-earth  environment.  Only  the  trapped  proton 
environment  lias  not  been  Included,  since  It  has  already  been  described  by 
the  AP-8  model  of  Sawyer  arid  Vette  (1976). 

Sections  2.0,  3.0  and  4.0  present  the  aata  base  for  particles  In  the 
1 nterplanetary  medium  and  describe  how  this  environment  has  been  modeled. 

Section  5.0  discusses  the  geomagnetic  cutoff  and  describes, 
qualitatively,  how  It  modifies  the  particle  spectra  from  the  Interplanetary 
medium.  The  second  report  of  this  study  will  describe  an  accurate  method 
for  modulating  the  Interplanetary  spectra  to  obtain  the  orbit-averaged 
spectra  Incident  on  any  spacecraft  In  any  orbit. 

The  composition  of  particles  trapped  in  the  earth's  magnetosphere  Is 
described  In  Section  6.0.  The  heavy  Ion  composition  of  trapped  radiation  at 
energies  above  10  MeV/u  is  the  least  well  known  part  of  the  particle 
environment.  The  few  measurements  that  exist  show  heavy  Ion  fluxes  higher 
than  those  In  the  Interplanetary  medium. 

Section  7.0  discusses,  In  a  qualitative  way,  how  shielding  alters  the 
particle  spectra.  Cosmic  ray  transport  theory  In  condensed  media  will  be 
the  subject  of  a  thlra  report.  This  section  also  reviews  the  work  that  has 
been  publlsljecj  t.o  date  on  soft  errors  and  gives  a  general  discussion  of  the 
environment  add  Its  effects  on  electronics  In  various  orbits. 

The  status  of  the  data  base  for  this  first  part  of  the  study  Is  reviewed 
In  Section  8.(J  and  recommendations  are  made  for  additional  work  that  would 
allow  the  particle  environment  model  to  be  Improved. 
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Fiy.  1.3  —  The  universal  abundance  of  the  chemical  elements  in  nature  relative 
to  silicon  *  106.  These  results  are  obtained  from  studies  of  meteorite*,  our  sun 
and  other  stars  (taken  from  Cameron,  1980). 


2.0  The  Galactic  Cosmic  Ray  Model 


Most  of  the  energetic  charged  particles  usually  found  In  the  vicinity  of 
earth  are  cosmic  rays,  particles  which  come  from  outside  our  solar  system. 
The  sources  of  these  cosmic  rays  are  as  yet  unknown.  The  existing  evidence 
suggests  that,  except  for  the  highest  energies,  these  particles  come  from 
sources  within  our  galaxy.  Cosmic  rays  arriving  In  our  solar  system  consist 
of  the  nuclei  of  J1  the  elements  In  the  periodic  table  and  electrons. 

2.1  The  Nucleonic  Component  of  Cosmic  Rays 

By  studying  the  differences  between  the  chemical  composition  of 
nucleonic  cosmic  rays  and  material  in  our  universe  generally,  we  have  been 
able  to  learn  some  things  about  the  cosmic  ray  population  In  our  galaxy.  As 
cosmic  rays  travel  through  the  galaxy,  they  occasionally  collide  with  nuclei 
of  Interstellar  gas.  The  resulting  nuclear  reactions  modify  the  composition 
of  cosmic  rays.  A  detailed  recent  estimate  of  these  modifications  has  been 
given  by  Sllberberg  et  al  (1S76).  These  authors  have  found  that,  by 
assuming  cosmic  rays  traverse  ~  5.5  g/cmz  of  Interstellar  gas  on  the 
average,  they  can  account  for  almost  all  the  differences  In  chemical 
composition.  These  results  are  further  supported  by  the  measured  cosmic  ray 
abundances  of  electron-capture  Isotopes  such  as  'Be  that  could  only  have 
been  produced  In  collisions  with  Interstellar  gas  (see  for  example, 
Wledenbeck  and  Greiner.  1580).  By  measuring  the  cosmic  ray  abundance  of  the 
radioactive  Isotope  10Be  (T1//z  ■  1.6  x  106  years),  Wledenbeck  and 
Crelner  (1980)  have  shown  that  cosmic  rays  reaching  earth  have  wandered 
about  In  our  galaxy  for  -  8  x  106  years,  on  the  average.  Their 
measurements  are  consistent  with  the  results  of  a  number  of  earlier 
Investigators. 

These  results  and  others  have  led  to  a  model  for  cosmic  ray  confinement 
In  the  galaxy.  The  standard  model  assumes  that  the  galaxy  Is  uni  formally 
populated  with  cosmic  ray  sources.  These  sources  emit  cosmic  rays  Into  the 
galaxy  where  they  diffuse  through  the  random  magnetic  fields  of  the  galaxy, 
but  are  contained,  with  some  leakage  at  the  galactic  boundary. 

In  the  context  of  the  standard  model,  Figure  2.1  shows  how  the  cosmic 
ray  composition  Is  tran-  formed  by  fragmentation  as  cosmic  rays  wander 
through  the  galaxy  on  their  way  to  earth,  Adams,  et  al.  (1980a).  The 
abundances  at  earth  are  plotted  on  a  scale  relative  to  arriving  carbon  - 
100.  The  abundances  are  broken  down  according  to  the  fraction  that  have 
survived  collisions  (open  bars)  with  Interstellar  gas  to  reach  us  and  those 
that  were  produced  by  collisions  of  heavier  cosmic  rays  with  Interstellar 
gas  (filled  bars).  Also  shown  are  the  inferred  source  abundances  (dashed 
bars).  It  should  be  noted  that  about  half  of  the  heavy  (Z  >  6)  cosmic  ray 
nuclei  have  collided  with  Interstellar  gas  nuclei. 

The  most  abundant  element  In  cosmic  rays  Is  hydrogen.  Figure  2.2  shows 
the  differential  energy  spectrum  of  hydrogen  (for  the  most  part  protons). 

The  data  shown  In  this  figure  are  only  the  most  recent,  measurements  of 
cosmic  ray  protons.  They  are  consistent  with  the  much  larger  number  of 
measurements  carried  out  In  the  50 1 s  and  60 1 s .  We  have  selected  the  data 

presented  below  10  KeV  to  show  only  those  measurements  made  during 
periods  of  maximum  and  minimum  solar  activity.  The  smooth  solid  curve  Is  an 
analytic  function  fit  to  the  data. 
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At  the  highest  energies,  the  proton  spectrum  has  the  mathematical  form 
of  a  power  law,  i.e.  dJ/dE  ~  £-<*,  with  a  spectral  index,  a  ~  2,75. 

Power  law  spectra  can  be  produced  by  particle  acceleration  in  random  moving 
magnetic  fields  as  shown  by  Fermi  (194??).  The  conditions  for  Fermi 
acceleration  occur  In  a  variety  of  astrophyslcal  settings.  We  believe  that 
outside  the  solar  system,  cosmic  rays  obey  a  power  law  to  much  lower 
energies  than  shown  in  Figure  2,2.  The  deviation  from  a  power  law  below 
5000  MeV/u  In  Figure  2.2  Is  largely  due  to  solar  modulation.  The  power-law 
fit  Is  better  In  the  case  of  a  rigidity  spectrum;  some  of  the  deviation  Is 
due  to  the  transformation  from  rigidity  to  kinetic  energy.  To  reach  the 
vicinity  of  earth,  cosmic  rays  must  "swim"  upstream  In  the  solar  wind.  The 
process  of  diffusion  Inward  against  the  outward-flowing  solar  wind  (see  Webb 
and  Gleason,  1SCC;  Jokipii,  et  al.  1977;  Fisk,  1976;  Jokipii,  1971)  reduces 
the  energy  of  the  cosmic  rays  an  average  ~  300-400  MeV/u.  It  also 
attenuates  the  flux  arriving  near  earth  In  an  energy  dependent  way.  The 
amount  of  solar  modulation  depends  on  the  general  level  of  solar  activity. 
When  the  sun  Is  quiet  and  especially  during  the  minimum  of  the  11-year  solar 
activity  cycle,  cosmic  rays  have  the  easiest  access  to  the  earth's  orbit. 
These  periods  account  for  the  upper  branch  of  the  spectrum  In  Figure  2.2. 

The  lower  branch  corresponds  to  a  quiet  (no  flares)  period  during  the 
maximum  of  the  11-year  cycle.  Solar  modulation  Is  a  complex  subject  and 
still  an  area  of  active  research.  We  will  describe  our  model  for  dealing 
with  It  In  a  later  section. 

At  the  lowest  energies  shown  in  Figure  2.2,  the  cosmic  ray  flux  varies 
considerably  even  when  no  large  solar  flare  is  In  progress.  These 
variations  take  the  form  of  short  term  increases  above  a  lower  limit  that 
varies  slowly  with  the  solar  cycle.  These  Increases  are  due  to  small  solar 
flares,  flares  poorly  connected  to  the  earth  (I.e,  on  the  backside  of  the 
sun)  and  particles  accelerated  by  the  solar  wind  In  co-rotating  Interaction 
regions  (CIR)  In  the  Interplanetary  medium  (to  be  discussed  In  section  3.0). 

These  variations  have  been  observed  on  IMP -8,  an  interplanetary  probe 
orbiting  the  earth  at  24  to  26  earth  radii.  Figure  2.3  shows  six-hour 
averages  of  the  proton  flux  observed  on  IMP-8  (Pyle,  1981)  as  a  function  of 
flux  level.  These  protons  had  energies  between  11.24  MeV  and  29.75  MeV  and 
the  data  span  the  period  from  Oct  30,  1973  to  July  2,  1980.  The  most  common 
flux  level  measured  was  In  the  range  of  the  galactic  cosmic  ray  background 
(OCR)  and  corresponds  to  the  range  between  the  solar  minimum  and  solar 
maximum  spectra  In  Figure  2.2.  The  tall-off  In  measurements  below  this  flux 
level  Is  due  to  temporary  Increases  In  solar  modulation  called  Forbush 
Decreases  (Forbush,  1938).  Above  the  flux-level  of  galactic  cosmic  rays, 
there  Is  a  long  tall  extending  up  for  many  orders  of  magnitude.  The 
smallest  of  these  Increases  Is  due  to  the  addition  of  protons  from 
co-rotating  Interaction  regions,  Fan,  et  al .  (1965)  (also  discussed  In 
section  3.1)  or  small  solar  flares.  Flux  levels  observed  between  3  and 
60,000  protons/m2ster  sec  MeV/u  are  due  to  medium-size  flares  or  larger 
ones  that  were  poorly-connected  to  IMP-8  by  the  Interplanetary  magnetic 
fields.  The  flux  levels  above  this  range  are  due  to  large  flares  which  are 
treated  separately  In  section  4,0.  Also  shown  in  Figure  2.3  is  a  90  per 
cent  confidence  level,  that  Is  a  flux  level  which  was  exceeded  in  only  10 
per  cent  of  the  six-hour  Intervals, 

Figure  2.2  shows  a  worst  case  proton  spectrum  (with  90  per  cent 
confidence),  based  on  four  energy  intervals  spanning  the  range  11.24  MeV  _<  E 
<  94,78  MeV. 
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2.2  The  Relative  Abundances  of  Cosmic  Rays 

It  would  be  convenient  to  simply  scale  the  hydrogen  spectrum.  Figure 
2.2,  according  to  the  ratio  of  hydrogen  with  respect  to  the  other  elements. 
Unfortunately  the  ratio  of  hydrogen  to  the  other  elements  depends  on 
particle  energy  and  the  level  of  solar  modulation.  Basically,  this  is 
because  the  charge  to  mass  ratio  for  hydrogen  Is  ~  1  while  it  Is  -  0.5 
for  the  other  elements.  This  leads  to  different  responses  of  the  spectra  to 
magnetic  rigidity  dependent  and  velocity  dependent  phenomena.  It  Is  better 
to  treat  hydrogen  as  a  special  case  and  proceed  to  helium. 

The  helium  differential  energy  spectrum  is  shown  in  Figure  2.4.  The 
data  points  shown  In  this  spectrum  are  only  a  representative  sample  of  the 
data  we  examined.  The  density  of  points  plotted  precluded  the 
Identification  of  each  data  point  with  Its  author.  To  avoid  cluttering  the 
figure  we  have  shown  error  bars  on  only  a  sampling  of  the  data  points.  The 
data  shown  come  from  measurements  made  throughout  the  solar  cycle,  though  we 
Included  as  many  data  points  as  possible  near  solar  maximum  and  minimum. 

The  helium  data  we  have  used  In  the  figure  came  from  Ryan  et  al.  (1972), 

Smith  et  al .  (1973),  Verma  et  al .  (1972),  Anand  et  al.  (1968),  Ormes  and 
Webber  (1965),  Von  Rosenvlnge  et  al.  (1969),  Webber  et  al.  (1973a),  Fan  et 
al .  (1965),  Balasubrahmanyan  et  al.  (1965),  Freler  and  Waddlngton  (1965), 
Hofmann  and  Winckler  (1966),  Cleghorn  et  al.  (1971),  Leech  and  O'Gallagher 
(1978),  Webber  and  Leznlak  (1973),  Bhatla  et  al.  (1977),  Rygg  and  Earl 
(1971),  Webber  and  Ormes  (1968),  Badhwar  et  al.  (1969),  Ormes  and  Webber 
(1968),  Balasubrahmanyan  et  al  (1967),  Mason  (1972)  and  Garcla-Munoz  et  al. 
(1975),  though  the  data  of  other  authors  was  also  consulted.  The  smooth 
solid  curve  is  from  an  analytic  function  we  have  fit  through  the  data  points. 

As  In  the  case  of  the  proton  spectra  discussed  earlier,  the  helium  flux 
levels  at  low  energies  are  sometimes  measured  to  be  considerably  different 
from  those  predicted  by  the  analytic  spectral  functions  In  Figure  2.4. 

Figure  2.5  shows  6-hour  averages  of  the  helium  flux  measured  on  IMP-8, 

(Pyle,  1981).  The  helium  nuclei  had  energies  of  10.9  MeV/u  <  E  <_  25.36 
MeV/u,  and  the  data  set  spans  the  same  period  as  the  proton  Bata  shown 
earl ler. 

During  most  of  the  period  covered  by  these  observations,  the  low  energy 
helium  spectrum  was  dominated  by  the  addition  of  anomalous  component  (to  be 
discussed  later  In  section  3.2).  This  accounts  for  the  location  of  the  most 
common  flux  level  measured  in  this  period.  Lower  flux  levels  were  measured 
after  the  spring  of  1978  when  the  anomalous  component  no  longer  contributed 
to  the  flux  near  earth.  These  two  conditions  are  smeared  together  by 
Forbush  decreases. 

As  in  the  case  of  protons  the  enhanced  flux  levels  are  the  result  of 
particles  accelerated  In  co-rotatlng  interaction  regions  (CIR's)  and  flares 
of  varying  sizes. 

Figure  2.4  shows  a  worst  case  spectrum  (with  a  90  per  cent  confidence 
level).  This  spectrum  Is  chosen  so  that  fluxes  above  this  level  are 
observed  only  10  per  cent  of  the  time.  These  data  are  based  on  four  energy 
Intervals  between  10.9  MeV/u  and  94.81  MeV/u. 
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Comparing  Figures  2.2  and  2.4  we  see  that  the  cosmic-ray  He  abundance  Is 
'  15  per  cent  of  the  H  abundance  In  the  energy  range  20C-700  MoV/u,  and 
**  5  per  cent  above  104  fieV/u.  Helium  is  the  best  element  to  choose  for 
measuring  the  differential  energy  spectrum  because  It  Is  distinct  from  all 
the  singly  charged  particles  (l.e.  protons,  electrons,  muons,  and  plons  all 
have  one  charge);  it  is  plentiful;  and  it  has  a  charge  to  mass  ratio  similar 
to  the  heavier  elements. 

As  mentioned  in  Section  2.1,  cosmic  rays  spend  ~  107  years  diffusing 
around  the  galaxy  and  being  broken  up  in  collisions  with  the  Interstellar 
gas.  Not  surprisingly,  this  diffusion  process  Is  energy  dependent  and  the 
higher  energy  cosmic  rays  have  not  travelled  as  far,  or  as  long  as  the  lower 
energy  ones.  This  means  that,  at  higher  energies,  there  will  be  fewer 
arriving  secondaries  and  more  surviving  primordial  cosmic  rays. 

Cosmic  ray  helium  is  mostly  surviving  primordial  material  In  the  context 
of  Figure  2.1;  only  ~  10  per  cent  of  He  is  secondary.  This  places  it  in 
the  same  class  with  hydrogen,  carbon,  oxygen,  neon,  magnesium,  silicon, 
sulfur,  calcium  and  iron;  a  list  which  includes  the  most  abundant  nuclei. 

We  would  expect,  as  Figure  2.6  shows,  that  the  helium  to  carbon  plus  oxygen 
ratio  is  nearly  independent  of  energy  at  21  +  2  for  the  1-5  GeV/u  range.* 
Figure  2.7  (from  Caldwell,  1977)  shows,  base?  on  less  data,  that  the  ratio 
of  (neon  +•  magnesium  +  si  1 1con)/hel lum  does  not  vary  much  with  energy. 

Figure  2.8  however  shows  that  the  Fe/He  ratio  is  energy  dependent.  To  some 
extent,  this  merely  reflects  the  relatively  larger  fraction  of  surviving  Fe 
at  high  energies.  In  this  way,  the  ratio  can  increase  by  a  factor  of  ~  3 
as  can  be  inferred  from  Figure  2.1.  The  Fe/He  ratio  could  increase  even 
more,  if  the  source  spectra  of  Fe  and  He  differ  as  well.  It  appears  that  Fe 
will  have  to  be  treated  separately  from  helium.  Figure  2.9  shows  the 
differential  energy  spectrum  of  Fe.  The  data  base  for  the  Fe  spectrum  is 
rather  limited.  Figure  2.9  shows  all  the  published  data  for  iron  from  10 
MeV/u  to  100  MeV/u.  Between  100  MeV/u  and  103  MeV/u  there  is  an  adequate 
set  of  measurements  during  solar  minimum  conditions,  but  there  are  no 
published  measurements  during  solar  maximum  (experiments  are  in  progress  at 
NRL  and  elsewhere  to  obtain  these  data).  For  the  present,  we  have  used  the 
general  shapes  of  the  solar  maximum  and  minimum  helium  spectra  as  a  guide  to 
obtain  the  smooth  solid  curve  shown  in  Figure  2.9.  By  analogy  with  the  flux 
enhancements  found  for  helium,  we  have  suggested  a  worst  case  spectrum  for 
iron  shown  as  a  dashed  line. 

The  differential  energy  spectra  of  all  the  elements  between  helium  and 
nickel  will  be  obtained  by  scaling  the  helium  or  iron  spectra.  Figure  2.10 
shows  the  data  on  elemental  composition  of  lithium  through  sulfur, 
normalized  to  helium  *  1000.  The  data  in  Figures  2.6  and  2.7  together  with 


*If  the  He,  C  and  0  source  spectra  are  identical,  this  ratio  is  -  15 
for  a  path  length,  X  <  1  g/cm2,  i.e.  for  E  >  50  GeV/nucl .  and  may  go  to 
-  23  for  X  "  8  g/cm2,  which  is  plausible  at  energies  of  200  to  600 
MeV/u.  We  have  adopted  a  ratio  of  21,  for  the  complete  integral  energy 

spectrum. 
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that  of  Juliusson  (1974),  Lezniak  and  Webber  (1978),  Orth  et  al  (1978)  and 
Caldwell  (1977)  show  that  the  ratios  of  C,  0,  Me,  Mg,  and  SI  to  He  are 
approximately  Independent  of  particle  energy.  This  is  to  be  expected  since 
as  Figure  2.1  shows  these  nuclei  are  principally  surviving  primordial 
material.  Table  2.1  shows  the  relative  abundances  we  have  adopted  for  these 
elements  as  well  as  sulfur. 


TABLE  2.1  The  Elemental  ratios  for  Elements  having  Helium-like 
and  Iron-like  Spectra  Respectively 


Element 

Ratio 

to  He 

Element 

Ratio 

to  Fc 

C 

2.5  x 

iO-2 

Ca 

2.3  x 

10"1 

0 

2.3  x 

10-2 

Co 

6.0  x 

10~3 

F 

4.1  x 

M 

o 

1 

-c* 

Ni 

4.8  x 

10'2 

Ne 

3.5  x 

10"3 

Na 

7.0  x 

10“4 

Mg 

4.7  x 

10'3 

Al 

8.3  x 

IO"4 

51 

3.5  x 

10“3 

P 

2.0  x 

10"4 

S 

7.4  x 

10'4 

Figure  2,11  shows  the  energy  dependence  of  the  ratio  of 
(Li  +8e  +  B)/He.  Since  all  three  elements  Li,  Be  and  B  are  entirely 
secondary,  we  believe  that  each  of  them  displays  this  energy  dependence. 
Table  2.2  shows  the  ratios  Li/(L1  +  Be  +  B),  Be/(L1  +  Be  +  B)  and  B/(Li  +  Be 
+  B).  Using  these  ratios,  we  can  scale  Figure  2.11  to  obtain  the  ratios 
Ll/He,  Be/He  and  B/He  as  a  function  of  energy.  Yhe  differential  energy 
spectra  of  these  elements  can  then  be  obtained  from  Figure  2.4. 

The  ratio  N/He  Is  shown  In  Figure  2.12.  It  is  also  clearly  energy 
dependent,  but  In  a  different  way.  Figure  2.12  car  be  used  to  scale  Figure 
2.4  to  obtain  the  nitrogen  differential  energy  spectrum. 

From  Figure  2.1,  we  can  guess  that  F,  Na,  Al  and  P  will  also  have  energy 
dependent  ratios  relative  to  He.  The  available  experimental  data  (see  Orth 
et  a1.,  1978;  Juliusson,  1974;  and  ^ezniak  and  Webber,  1978}  are  not  of 
sufficient  accu.acy  to  define  this  energy  dependence,  so  we  will  use 
constant  ratios.  The  adopted  values  are  shown  in  Table  2.1 

Figure  2.13  shows  the  ratios  of  the  elements  17  £  Z  <  25  to  Fe  as  a 
function  of  energy.  While  Tfil s  ratio  is  energy  dependent,  it's  not  clear 
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that  all  the  elements  In  the  numerator  display  this  dependence.  Figure  2.1 
shows  that  calcium  is  mostly  primordial  material,  we  would  therefore  expect 
Ca/fe  to  be  Independent  of  energy.  The  abundances  of  the  elements  Cl 
through  hi  are  shown  in  Figure  2.14  normalized  to  Fe  =  100.  The  adopted 
value  for  Ca/Fe  in  shown  In  Table  2.1 

The  adopted  ratios,  at  low  energies,  of  the  other  elements  in  the  17  1 
1  1  25  range  to  the  sum  of  the  elements  in  that  range  are  shown  In  Table 
2.2.  These  ratios  are  used  to  scale  the  energy  dependent  ratio  In  Figure 
2.13  so  as  to  obtain  the  ratios  Cl /Fe ,  etc.  which  in  turn  are  used  to  scale 
the  Fe  spectrum,  Figure  2.9,  to  the  spectra  of  these  elements. 


TABLE  2.2  The  Elemental  Ratios  Required  to  obtain  the 
Individual  Elemental  Spectra  from  Figures  2.11  and  2.13 
Combined  with  Figures  2.4  and  2.9  respectively. 


Ratio 

Relative 

Abundances 

L1/(L1  +  Be  +  B) 

0.33 

Be/(L1  +  Be  +■  B) 

0.175 

B/(L1  +  Be  +  B) 

0.5C 

01/(17  <  2  <  25) 

0.07 

Ar/ (17  <l<iS) 

0.13 

K/U7  <  Z  <  25) 

0.09 

Sc/ (17  <  Z  <  25) 

0.05 

T1/(17  <Z  <  25) 

0.14 

V/(17  <  Z  <  25) 

0.07 

Cr/( 17  <  Z  <  26) 

0.14 

Mn/{ 17  <2  <  25) 

0.10 

2.3  Nuclei  Heavier  than  Nukel 


The  galactic  cosmic  rays  consist  of  every  element  in  the  periodic 
table.  So  far  we  have  only  dealt  with  the  first  28,  which  are  the 
most  abundant.  The  abundances  of  the  remaining  elements  relative  to 
106  Fe  are  shown  In  Table  2.3,  (Adams,  et  a  1 . ,  1980b). 
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TABLE  2.3  Abundances  of  Trans-Iron  Nuclei  in  Galactic  Cosmic  Rays 


Atomic  Number 
26 

29  <  Z  <  34 
35  <  Z  <  39 
Z  >  40 


Relative  Abundance 
1C6 

1.7  x  103 
1.7  x  102 
8  x  101 


For  the  purposes  of  this  particle  environment  model,  these  very  rare,  but 
very  damaging  nuclei  will  be  ignored.  It  should  be  noted,  however,  that 
should  a  microelectronic  component  be  struck  by  one  of  these  rare  nuclei,  an 
enormous  amount  of  charge  would  be  liberated,  leading  to  a  soft  error  even  In 
devices  commonly  thought  to  be  Insensitive  to  this  effect! 

2.4  Cosmic  Ray  Electrons 

There  appear  to  be  two  plausible  methods  by  which  electrons  can  produce 
soft  errors.  The  first  Is  by  directly  depositing  enough  energy  In  the 
critical  volume  to  produce  the  required  critical  charge.  The  second  Is  by 
producing  bremstrahlung  photons  that,  in  turn,  undergo  photo-nuclear 
Interactions  with  the  silicon  In  the  device. 

We  will  consider  the  direct  method  first.  Electrons  deposit  energy  most 
densely  near  the  end  of  their  range.  Because  of  their  low  rest  mass, 
electrons  undergo  large  angle  scattering  before  their  stopping  power  has 
risen  much  above  Its  minimum  value.  This  causes  the  practical  range 
(displacement  distance)  of  a  stopping  electron  to  be  much  shorter  than  Its 
path  length  with  the  result  that  the  electron  deposits  all  Its  energy  In  a 
relatively  small  volume.  The  practical  range  of  an  electron  in  aluminum  Is 
given  by: 

r  =  5.37  x  10-lECl-O. 9815/(1  -  3.123E)]  g/cm2  (2.1) 

where  E  Is  in  MeV  (see  Kobetlch  and  Katz,  1968).  Without  Introducing  much 
error  we  may  use  this  equation  for  silicon  and  compute  the  electron  energy 
corresponding  to  any  practical  range.  If  this  practical  range  Is  taken  to  be 
the  diameter  of  a  collection  volume,  then  the  corresponding  energy  Is  roughly 
the  energy  one  might  expect  an  electron  to  deposit  In  that  volume.  Using  3.6 
ev  per  electron-hole  pair,  we  can  estimate  the  charge,  Q,  that  the  electron 
produces. 

Figure  2.15  shows  Q  In  electron-hole  pairs  as  a  function  of  the  mean 
device  diameter.  This  figure  suggests  that  devices  such  as  the  256K  CCD 
described  by  Ziegler  and  Lanford  (1979)  vdl1  have  soft  errors  due  to  stopping 
electrons.  It  should  be  noted  that  these  n^ed  not  be  cosmic  ray  electrons; 
trapped  electrons,  air  shower  electrons  and  electrons  from  terrestrial  r-ray 
Interactions  would  be  equally  effective! 

While  electrons  seem  to  be  capable  of  producing  errors  directly  In 
devices  sensitive  to  <  104  electron-hole  pairs,  they  are  unable  to  produce 
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errors  when  >  10$  electron-hole  pairs  are  required.  Devices  being 
considered  for  satellite  applications  are  much  less  sensitive  than  256K  CCD's 
and  cannot  be  directly  upset  by  stopping  electrons. 

The  second  way  In  which  electrons  can  cause  errors  Is  effective  for  less 
sensitive  devices.  Electrons  must  produce  photons  that,  In  turn,  undergo 
SI ( y ,  n),  SI ( y .  p),  or  SI  ( y , a )  reactions  In  the  devices.  Because  of  the 
thresholds  for  these  reactions,  electrons  with  energies  below  20  MeV  will  not 
cause  these  reactions.  Webber  (19735  has  reviewed  the  cosmic  ray  electron 
differential  energy  spectra.  He  shows  that  the  electron^flux  Is  comparable 
to  the  proton  flux  at  10  MeV,  but  falls  rapidly  to  -  10"2  of  the  proton 
flux  at  100  MeV.  Clearly,  low  energy  protons  produced  by  electrons  will 
always  be  out-numbered  by  cosmic  ray  protons.  As  was  shown  In  section  2.2, 
the  alpha  flux  is  ~  15  per  cent  of  the  proton  flux,  so  electron-produced 
alpha  particles  will  always  be  overwhelmed  by  cosmic  ray  alphas. 

In  general,  we  conclude  that  low  energy  electrons  (<  20  MeV)  will  not 
cause  errors  In  the  relatively  Insensitive  components  considered  for 
satellite  applications.  Higher  energy  electrons  can  cause  errors  by  the 
three  stage  process  described  above,  but  this  process  will  be  Important  only 
If  the  electron  flux  Is  enormously  larger  than  the  elemental  flux. 

2.5  Solar  Modulation 


As  can  be  seen  In  Figures  2.2,  2.4,  and  2.9,  the  differential  energy 
spectra  are  spread  between  two  extremes  below  -  103  MeV.  This  Is  due  to 
solar  modulation  of  the  differential  energy  spectra  Incident  on  the  solar 
cavity  and  depends  on  the  level  of  solar  activity. 

Figure  2.16  shows  the  annual  average  cosmic  ray  flux  for  the  past  four 
decades,  measured  for  most  of  that  period  by  the  neutron  monitor  at  Deep 
River  (Rao,  1972,  and  Ahluwalla,  1979).  This  monitor  detects  hadrons, 
primarily  neutrons,  which  are  secondary  products  of  cosmic  rays  incident  on 
the  atmosphere.  In  this  way  It  measures  the  cosmic  ray  flux  at  earth 
continuously.  The  valleys  In  1947,  1958  and  1969  correspond  to  maxima  In 
solar  activity.  The  detailed  shape  of  the  curve  over  several  solar  cycles  Is 
quite  variable,  though  crudely  sinusoidal. 

To  estimate  the  low  energy  spectra  at  any  time  in  the  past,  It  Is  best 
to  peg  the  modulation  level  by  the  measured  intensity  In  experiments  carried 
out  at  that  time  or,  at  best  another  time  when  the  solar  neutron  monitor 
levels  were  similar.  The  solar  modulation  level  In  the  near  future  may  also 
be  predicted  by  extrapolating  the  present  solar  neutron  monitor  level,  using 
a  sine  curve  with  the  same  period  as  that  shown  in  Figure  2.16.  This  method 
Is  probably  only  reliable  for  predictions  less  than  on*')  year  Into  the  future. 

In  modeling  the  spectra  of  cosmic  rays  for  satellite  planning,  we  must 
be  able  to  predict  the  level  of  solar  modulation  years  Into  the  future.  It 
seems  that  a  simple  sine  function: 

M  *  A  sin  W(t-t0 )+B  (2.2) 

Is  the  best  choice.  The  function,  eq.  (2.2)  with  W  ■  2ir/10.9  years  ■ 
0.576/year  and  t0  *  1950.06  Is  shown  as  the  smooth  curve  In  Figure  2.16. 
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The  values  of  A  and  B  were  chosen  to  best  fit  the  data  In  Figure  2.16.  As  we 
can  see,  the  cosmic  ray  rlux  Is  quite  varlaole  from  one  solar  cycle  to  the 
next,  and  only  crudely  predicted  by  eq.  2.2.  The  Deep  River  neutron  monitor 
responds  mostly  to  very  energetic  cosmic  rays,  so  the  amplitude  of  the  solar 
cycle  variation  Is  much  greater  at  lower  energies  (see  Figure  2.2,  2.4  and 
2.9),  and  probably  less  predictable.  We  feel  that  our  present  Inability  to 
predict  the  level  of  solar  modulation  In  the  future  Is  the  principal  source 
of  uncertainty  In  the  estimates  of  future  cosmic  ray  flux  levels  below 
-  1000  MeV/u. 

2.6  The  Analytic  Model  for  Galactic  Cosmic  Rays 

In  the  preceding  sections  we  have  discussed  the  nature  of  the  cosmic  ray 
energy  spectra  and  chemical  composition.  In  this  section  we  present  a  simple 
analytic  recipe  that  may  be  used  to  estimate  the  differential  energy  spectrum 
of  each  of  the  first  28  elements  In  cosmic  rays  between  10  and  105  MeV/u. 

The  solid  curves  In  Figures  2.2,  2.4  and  2.9  are  chosen  to  form  an 
envelope  around  the  data  at  low  energies  and  blend  to  a  single  curve  at  high 
energies.  The  analytic  form  of  these  Is: 

F(E,t)  -  A(E ) sin  w(t-t0)  +  8(E)  In  partlcles/m^ster.sec.MeV/u  (2.3) 
where  W  ■  0.576/year,  t0  ■  1950.6, 

A(E)  ■  0.5[fm1n  (E)  -  f„,ax  (E)], 

B(E)  •  O.BCfmin  (E)  +  fniax(E)] 

The  spectral  shapes  fmaX  and  fmin  are  both  obtained  from  the 
equation: 

f(E)  -  ion  (E/E0)a  (2.4) 

where 

a  ■  a0  |  1  -  exp[-Xi(logioE)b3}  (2.5) 

and 

m  ■  Ci  expC-X2(lugioE)^3“C2  (2.6) 

The  constants  a0,  E0l  b,  Xj.,  X2,  Ci  and  C2  are  given  In  Table  2.4 

for  the  solar  maximum  and  solar  minimum  cases  or  the  proton,  helium  and  Iron 

spectra. 
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TABLE  2.4  Parameter  Values  used  In  Eq.  (2.4)  to  Reproduce  the  solar 
maximum  and  solar  minimum  envelopes  (solid  curves)  shown  in 
Figures  2.2,  2.4  and  2.9  for  hydrogen,  helium  and  iron  respectively 


Element 

Solar 

Activity 

ao 

Eo 

b 

*1 

*2 

Cl 

C2 

hydrogen 

min 

-2.2 

117500 

2.75 

.117 

.80 

6.52 

4.0 

hydrogen 

max 

-2.2 

117500 

2.75 

.079 

,80 

6.52 

4.0 

helium 

min 

-2.2S 

79400 

2.30 

.22 

.83 

5.0 

5.0 

helium 

mux 

-2.25 

79400 

2.30 

.155 

.83 

5.0 

5.0 

Iron 

min 

-2.70 

110000 

2.30 

.140 

.65 

7.0 

8.0 

Iron 

max 

-2.70 

110000 

2.30 

.117 

.65 

7.0 

8.0 

As  discussed  In  section  2.2,  flux  levels  well  above  the  cosmic  rty 
background  level  often  occur  at  energies  below  100  MeV/u.  We  recommend  that 
these  results  be  treated  as  follows:  1)  large  solar  flares  be  treated  as 
random  events  obeying  the  probability  distribution,  spectra  and  compositions 
discussed  In  Section  4.0;  2)  smaller  enhancements  be  handled  on  a  worst-case 
basis,  using  the  worst  case  spectra  (dashed  curves)  In  Figures  2.2,  2.4  and 
2. S’.  These  are  obtained  from  the  solar  minimum  cosmic  ray  spectra  discussed 
above.  For  protons, 

F worst  (E)  -  fm1n(E)Ll897e-E/9.66  +  1.64] 

For  helium  and  Iron  nuclei, 

E worst  (E)  -  fm1n(E)[28.4e-E/13.84  +  1.64] 

The  spectra  for  the  elements  C,  0,  F,  Ne,  Na,  Mg,  A1 ,  Si,  P  and  S  are 
all  obtained  by  scaling  Eq.  2.3  for  helium.  That  is,  just  compute  the  helium 
spectra  for  solar  maximum  and  solar  minimum  using  Eqs.  2.3,  2.4  and  Table 
2.4,  then  multiply  the  result  by  the  appropriate  entry  In  Table  2.1.  In  the 
same  manner,  the  spectra  of  the  elements  Ca,  Co  and  N1  are  all  obtained  from 
Eqs.  2.3  and  2.4  evaluated  for  Iron  and  multiplied  by  the  appropriate  entry 
from  Table  2.1. 


The  spectra  of  the  remaining  elements  are  more  complicated  to  obtain 
since  they  Involve  energy  dependent  charge  ratios.  Figure  2.11  shows  the 
ratio  of  (LI  +  Be  +  B)  to  He.  The  smooth  curve  In  the  figure  Is  the  ratio  we 
have  adopted.  Specifically  the  helium  spectrum  Is  modified  as  shown  below  to 
obtain  the  (LI  +  Be  +  B)  spectrum: 


F* 


/  0.0142  FHe,  E  <  6  x  103  MeV/u 


(2.7) 


I  0.67E-0.443  FHe,  E  >  6  x  103  MeV/u 

The  spectra  of  LI,  Be  and  B  are  obtained  by  multiplying  Eq.  (2.7)  by  the 
appropriate  ratios  In  Table  2.2 


1 


i 


t 

Jt 
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The  nitrogen  (N)  spectrum  must  also  be  obtained  by  an  energy  dependent 
modification  of  the  helium  spectrum,  Figure  2.12,  l.e. 

Fn  -{6.4  x  10“3  exp[-0.4(logioE-3.15)2]  +  5.6x10-3 

o  .  (2>6) 

expC-0.9{log3oE-0.8)23  f  FHe 


where  E  Is  In  MeV/u. 

The  spectra  of  the  elements  chlorine  (Cl),  argon  (Ar),  potassium  (K), 
scandium  (Sc),  titanium  (T1),  vanadium  (V).  chromium  (Cr)  and  manganese  (Mn) 
are  all  obtained  by  modifying  the  Iron  (Fe)  spectrum  with  the  function  Q(E)  to 
obtain  the  spectrum  of  these  combined  elements,  l.e. 

Q(E)  ■  16[1  •  exp(-.126E0.4)]E-.33 


^comb  ■  Q(E)  Flron^E) 

The  spectra  of  the  Individual  elements  are  obtained  by  multiplying  Fcomb  by 
the  appropriate  entry  from  Table  2.2. 

If  there  Is  an  Interest  In  the  spectra  of  the  elements  heavier  than 
nickel,  they  can  be  obtained  by  multiplying  FirQn,  Eqs.  2.3  and  2.4,  by  the 
appropriate  entry  In  Table  2.3. 


z 

Fig.  2.1  —  This  figure  shows  how  the  cosmic  ray  composition  is  transformed  by 
fragmentation  as  the  coBmic  rays  we  detect  wander  through  the  galaxy  on  their 
way  to  earth.  The  dashed  bars  show  the  source  composition.  The  adjacent  bars 
show  the  arriving  composition  with  the  filled  portion  of  the  bar  being  nuclei 
produced  by  fragmentation  and  the  open  portion  surviving  primordial  nuclei. 
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Fig.  2.2  —  The  differential  energy  spectrum  of  hydrogen  (mostly  protons).  The  data 
are  selected  to  show  the  solar  maximum  and  solar  minimum  forks.  The  smooth  curve 
is  an  analytic  function  contrived  to  fit  the  data.  The  duhed  curve  is  a  worst-case 
spectrum. 
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Fig.  2.3  •-  A  sample-size  distribution  of  six-hour  average  flux  measurements  of  protons 
(11.24  ivleV  <  E  <  29. 75  MeV),  The  measurements  were  made  by  the  University  of 
Chicago  Instrument  on  IMP-8  at  24  to  28  earth  radii  between  Oct.  30, 1973  and  July  2, 
1980  (Pyle,  1981),  Shown  among  the  bottom  are  the  sources  of  the  particles  that  con¬ 
tributed  to  the  flux  levels  above.  In  addition  to  solar  flare  protons,  those  sources  include 
galactic  cosmic  rays  (GCR)  and  protons  accelerated  in  co-rotating  interaction  regions 
(CTR). 
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Fig.  2 A  —  The  differential  energy  spectrum  of  helium  (mostly  alphas).  The  data  are  so 
extensive  that  it  was  not  possible  to  individually  attribute  the  data  points.  The  smooth 
curve  is  an  analytic  function  contrived  to  fit  the  data  for  solar  maximum  and  solar  min¬ 
imum  conditions.  The  dashed  curve  u  a  worst  case  spectrum. 


Fig.  2.5  —  A  sample-size  distribution  of  six-hour  average  flux  measurements  of  helium 
nuclei  10.9  MeV/u<  E  <  28.36  MeV/u),  The  measurements  were  made  by  the  Univer¬ 
sity  of  Chicago  instrument  on  IMP-8  at  24-26  earth  radii  between  Oct.  30,  1973  and 
July  2, 1980  (Pyle,  1981).  AC  in  this  figure  refers  to  the  anomalous  component,  refer 
to  Fig.  2.3  for  additional  details. 


ENERGY  (MeV'u) 

Fig,  2.6  —  The  He/(C+Oj  ratio  as  a  function  of  energy.  This  ratio  is  nearly 
energy  independent,  permitting  us  to  scale  the  He  spectrum  to  obtain  carbon 
and  oxygen  spectra.  Some  of  the  data  points  are  based  on  measured  ratios  like 
He/CNO,  He/O,  or  He/C.  These  were  corrected  to  obtain  He /(C+O)  ratio.  We 
adopt  the  value  21  ±  2  for  this  ratio  at  all  energies. 
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Fig.  2.8  —  The  Fe/He  ratio.  These  data  had  to  be  converted  from  Fe/(C+0),  and 
similar  ratios  to  Fe/He  using  the  He/(C+0)  ratio  from  Fig.  2.6.  ThiB  figure  dem¬ 
onstrates  that  the  He  and  Fe  spectra  do  not  have  the  same  shape. 
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Fig.  2.9  —  The  differential  energy  spectrum  for  Fe.  This  spectrum  has  not  been  defined 
well  at  low  energies  and  no  measurements  exist  during  solar  maximum.  As  a  result  the 
smooth  curves  were  defined  using  the  data  shewn  here  ao  well  as  the  He/Fe  data  shown 
in  Fig.  2,8  and  the  helium  spectrum,  Fig.  2.4.  The  smooth  curve  is  an  analytic  function, 
describing  the  iron  spectrum  at  the  extremes  of  solar  activity,  The  dashed  curve  is  a 
worst-case  spectrum. 
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Fig,  2,11  —  The  (Li  +  Be  +  B)/He  ratio;  the  data  shown  here  are,  for  the  rr-3t  part 
ratios  to  C  +  0.  These  were  corrected  to  He  using  the  results  of  Fig.  2.6.  Thy  solid 
curve  is  an  analytic  fit  through  the  data. 
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Fig.  2.12  The  N/He  ratio  it  ihown  at  a  function  of  energy.  The  data  pointi  were 
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FIs.  2.13  —  The  sub-iron  to  iron  ratio  is  shown  m  a  function  of  energy.  Some 
of  the  data  points  had  to  be  corrected  to  this  ratio  from  similar  measured  ra¬ 
tios.  The  Rmooth  curve  is  an  analytic  fit  through  the  data. 
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Fig.  2.14  -  The  abundances  of  the  element  chlorine  through  nickel  are  ihown  relative 
to  iron.  Theae  data  were  uaed  to  establish  the  charge  ratios  shown  in  Tables  2.1  and 
2.2  and  the  total  arriving  abundances  in  Fig.  2.1. 
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Fla.  2. IB  —  Thin  figure  shows  the  number  of  hole-electron  pairs  typically 
generated  by  an  electron  coming  to  rest  in  a  collection  volume  of  silicon 
as  a  function  of  the  me*)  diameter  of  that  volume, 
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Fig.  2.18  ~  The  ll-y-ar  solar  cycle  ii  uhown  u  measured  by  the  ground  level  cosmic 
ray  Intensity  at  the  Deep  River  neutron  monitor  (located  in  northern  Ontario,  Can¬ 
ada).  These  dnta  are  taken  from  Rao  (1972)  and  Ahluwolia  (1979).  The  imooth  curve 
is  a  sinusoidal  fit  to  tho  data.  The  large  depression  in  1958  refloats  the  extraordinary 
level  of  solar  activity  during  tho  19th  solar  cycle,  The  depressions  in  1961  and  1974 
could  be  a  22-year  repetitive  feature,  but  the  data  ore  not  sufficient  to  conclude  thu. 
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3.0  Particles  from  the  Interplanetary  Medium 

There  are  two  known  components  of  the  energetic  charged  particle 
environment  which  appear  to  originate  In  the  interplanetary  medium.  We  shall 
now  discuss  what  is  known  about  these  components  and  how  they  contribute  to 
the  particle  environment. 

3.1  Co-rotating  Events 

The  best  established  component  from  the  Interplanetary  medium  is  the 
co-rotating  particle  stream.  These  streams  are  correlated  with  high-speed 
solar-wlna  streams  and  Interplanetary  magnetic  field  structures  co-rotating 
with  the  sun.  The  particles  are  thought  to  be  selected  from  the  high  energy 
tail  of  the  solar  wind  and  accelerated  to  higher  energies  In  one  of  several 
ways  (see  the  review  by  Gloeckler,  1P7L). 

These  events  are  Infrequent  ant  produce  modest  Increases  in  the  particle 
flux  up  to  -  20  heV/u,  therefore,  they  affect  only  the  lowest  energies  of 
Interest  iri  this  study.  These  events  are  the  source  of  part  of  the 
fluctuations  in  the  low  energy  cosmic  ray  spectra  discussed  in  section  2.0. 

3.2  The  Anomalous  Component 

A  much  more  important  contribution  to  the  energy  spectra  comes  from  the 
anomalous  component.  This  component  was  discovered  by  Garcla-Munoz  et  al . 
(1972)  and  Independently  by  Hovestadt  et  al .  (1973)  and  McDonald  et  al. 
(1974).  Figure  2.1  (taken  from  Gloeckler,  1979)  shows  these  unusual  spectral 
features.  The  helium  spectrum,  Instead  of  dipping  to  a  minimum  at  about  10 
MeV/u  as  do  the  proton  and  carbon  spectra,  is  almost  flat  from  if  MeV/u  to  200 
MeV/u.  Notice,  that  tho  helium  flux  actually  exceeds  the  proton  flux  from 
4  MeV/u  to  30  MeV/u ! 

The  oxygen  spectrum,  while  following  the  carbon  spectrum  down  to  •  30 
MeV/u,  has  a  huge  peak  from  1  MeV/u  to  20  MeV/u.  Several  explanations  have 
been  offered  for  these  unusual  spectral  features.  The  most  widely  accepted 
theory,  due  to  Fisk  et  al .  (1974),  suggests  that  these  particles  come  from 
neutral  Interstellar  gas  that  can  freely  enter  the  heliosphere.  This  gas 
becomes  singly  ionized  as  it  approaches  the  sun.  Once  ionized  it  is 
accelerated  in  collision  regions  between  fast  and  slow  moving  streams  of 
solar  wind.  Because  of  the  very  good  vacuum  in  Interplanetary  space,  these 
particles  will  remain  singly  ionized  regardless  of  the  energy  they  acquire. 

Tho  Fisk  theory  predicts  that  only  atoms  with  first  ionization  potentials 
higher  than  hydrogen  will  display  anomalous  spectra  and  that  the  anomalous 
particles  will  be  singly  ionized.  The  first  prediction  has  largely  been 
borne  out  by  experiments  that  have  shown  anomalous  spectra  for  He,  N,  0  and 
Me,  but  not  for  H,  LI,  Be,  B,  C  and  F.  Tests  of  the  second  prediction,  that 
the  particles  are  singly  ionized,  have  so  fa«"  been  only  indirect  and 
inconclusive.  An  experiment  being  prepared  at  the  Naval  Research  Laboratory 
(Adams,  et  al.  19fl0b)  will  use  the  earth's  magnetic  field  to  test  this 
prediction. 

The  anomalous  component  Is  not  always  present  in  th"  vicinity  of  the 
earth.  It  appeared  between  1971  anu  1972  and  disappeared  again  with  the 
return  of  solar  maximum  in  1978.  It  is  unclear,  f>*oin  data  taken  during  the 
last  solar  minimum,  whether  the  anomalous  component  was  present  then. 
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One  theory  of  solar  modulation  (Joklpll  et  al.,  1977)  suggests  that  the 
anomalous  component  will  appear  near  earth  only  once  every  other  solar 
minimum,  i.e.  not  again  until  ~  1994.  It  remains  to  be  seen  whether  the 
anomalous  component  reappears  ~  1983  or  not  for  another  11  years. 

The  anomalous  component  Is  observed  even  more  strongly  In  the  outer  solar 
system  by  yoyager  and  Pioneer  spacecraft  and  remains  present  beyond  ~  10 
earth  radii  even  now.  Probably,  a  small  part  of  the  flux  observed  at  higher 
energies  near  the  earth  originates  In  the  anomalous  component  even  during 
solar  maximum. 

If  the  anomalous  component  Is  singly  fonlzed,  it  will  be  able  to 
penetrate  much  more  deeply  Into  the  earth's  magnetosphere  than  galactic 
cosmic  rays  at  the  same  energy.  It  could,  therefore,  make  a  much  more 
Important  contribution  Inside  the  magnetosphere  than  It  does  In  the 
Interplanetary  medium. 

From  Figure  3.1,  It  Is  clear  that  the  largest  contribution  to  the 
particle  spectra  of  Interest  here  Is  to  the  helium  spectrum.  The  enhance¬ 
ments  In  tne  N  and  0  spectra  are  at  energies  below  30  MeV/u  In  the 
Interplanetary  medium.  Secondly,  if  the  particles  are  singly  Ionized, 
geomagnetic  filtering  could  make  the  anomalous  component  more  Important  In 
the  magnetosphere.  Since  the  charge  state  of  these  particles  has  not  been 
established,  we  don't  know  how  the  earth's  magnetic  field  affects  them. 

Until  their  charge  state  has  been  established,  we  will  assume  that  they  are 
fully  Ionized  for  the  convenience  this  offers  In  treating  them. 

Our  lack  of  knowledge  of  the  charge  state  Introduces  a  second 
uncertainty.  Assuming  that  the  anomalous  component  Is  singly  Ionized,  Blake 
and  Frlesen  (1977)  have  suggested  that  anomalous  nuclei  entering  the 
atmosphere  might  be  stripped  In  the  lower  geocorona,  thus  becoming  stably 
trapped  for  periods  up  to  a  year  or  more.  This  could  add  heavy  Ions, 
unexpectedly,  to  the  trapped  radiation.  Also,  this  component  might  persist 
long  after  the  anomalous  component  became  undetectable  In  the  Interplanetary 
medium  In  1973.  Since  no  experimental  tests  of  this  theory  have  been 
performed,  we  will  assume  for  the  present  that  these  effects  do  not  occur. 

The  contribution  of  the  anomalous  component  to  the  helium  spectrum  can  be 
Included  by  ass1  ilng  a  constant  flux  extending  down  from  the  peak  flux  of 
the  spectrum  to  10  MeV/n.  For  the  oxygen  spectrum,  we  use  the  smooth  curve 
fit  through  the  anomalous  oxygen  peak  In  Figure  3.1.  The  equation  for  this 
curve  Is: 

f(E)  »  6xlO-?expC-(ln  E  -  1. 79)2/0. 70jpart1cles/m2ster. sec. MeV/u  (3.1) 

with  E  In  MeV/u. 

This  equation  should  be  used  In  the  energy  Interval  10  MeV/u  <  E  <  30 
MeV/u  to  replace  that  segment  of  the  cosmic  ray  oxygen  spectrum  (Tee  Section 
2.0).  This  will  Include  the  anomalous  component  In  the  cosmic  ray  oxygen 


spectrum.  In  a  similar  manner,  the  anomalous  component  can  be  added  to  the 
cosmic  ray  nitrogen  spectrum.  Using  the  N/0  ratio  reported  by  Klecker  et  al • 
(1977),  the  cosmic  ray  nitrogen  spectrum  Is  replaced  by: 

f(E)  *  l.S4xlO-2expL'-(ln  E-l. 79)2/0. 70]particles/m2ster. sec. MeV/u  (3.2) 

In  the  10  MeV/u  <  E  _<  30  MeV/u  energy  Interval. 

We  recommend  that,  helium,  nitrogen  and  oxygen  cosmic  ray  spectra  be 
altered  to  Include  the  anomalous  component  only  for  periods  of  solar 
minimum.  The  next  such  period  Is  1983-1989. 
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Fig.  3.1  —  The  quiet-time  spectra  of  hydrogen,  helium,  oxygen  and  carbon  in  the 
interplanetary  medium  during  the  last  solar  minimum,  1972-1978,  (taken  from 
Gloeckler,  1979).  It  is  instructive  to  compare  the  helium  spectrum  given  here 
with  that  in  Fig.  2.4.  The  solar  minimum  spectrum  (upper  solid  curve)  in  Fig.  2.4 
decreases  gradually  from  200  MeV/u  to  a  minimum  at  ~  15  MeV/u.  In  contrast, 
the  helium  spectrum  in  this  figure  remains  nearly  constant  below  200  MeV/u, 
rising  slightly  below  ~  50  MeV/u  and  even  exceeding  the  proton  flux  below  ~ 
30  MeV/u.  Tliia  difference  ia  due  to  the  additional  flux  provided  by  the  anomalous 
component  during  the  last  solar  minimum. 

| 


36 


4.0  Solar  FI  are  Parti c 1 e  s 


Solar  flares  are  sudden  outbursts  on  the  visible  surface  (photosphere)  of 
the  sun  which  release  huge  amounts  of  energy.  Most  of  this  energy  Is 
radiation  in  UV  and  X-rays.  A  part  of  this  energy,  mostly  from  hard  X-rays 
goes  Into  very  rapid  heating  of  the  solar  corona  above  the  flare.  This 
produces  large  currents  and  moving  magnetic  fields  In  the  corona  that 
accelerate  ambient  coronal  material  to  very  high  energies  quickly.  (For  a 
review  of  solar  flare  particle  acceleration,  see  Ramaty  et  al.  1980). 

Many  of  these  coronal  particles  escape  the  sun  and  spray  out  Into  the 
interplanetary  medium.  As  the  particles  move  Into  the  Interplanetary  medium 
they  tend  to  bt  guided  along  the  existing  spiral  magnetic  field  pattern  In 
the  ecliptic  plane.  As  a  result,  both  the  Intensity  and  the  spectrum 
observed  at  earth  depend  on  the  relative  positions  of  the  earth  and  the  flare 
on  the  sun.  For  example,  a  solar  wind  velocity  of  430  Km/sec  produces  a 
spiral  field  that  connects  the  earth  directly  to  points  on  a  solar  longitude 

line  -  54'  west  of  the  center  of  the  sun  as  viewed  from  earth.  For  flares 
at  other  solar  longitudes,  the  earth  will.  In  general,  receive  a  smaller  flux 
of  solar  particles;  the  flux  will  build  up  more  slowly;  and  It  may  contain 
fewer  high  energy  particles.  The  actual  degree  of  "well  connectedness" 
between  the  earth  ana  the  flare  site  depends  on  interplanetary  conditions  at 
the  time  of  the  flare  and  these  conditions  are  highly  variable  and 
unpredictable.  This  effect  may  lead  to  variations  as  large  as  100  In  the 
observed  flux  from  the  same  flare  at  different  points  around  the  earth's 
orbit  (see  Slmnett,  1976). 

4.1  The  Sizes  and  Frequencies  of  Flares 

Major  solar  flares  occur  at  random,  with  a  frequency  that  varies  from  one 
every  two  months  to  one  every  two  years.  The  particle  events  near  earth  that 
result  from  these  flares  last  from  two  hours  to  ten  days.  The  result  Is  that 
~  98  per  cent  of  the  time  the  particle  environment  In  the  Interplanetary 
medium  near  earth  Is  determined  by  galactic  cosmic  rays,  possibly  enhanced  at 
lew  energies  by  small  flares  etc.  (Section  2.0)  and  with  a  possible 
contribution  from  the  anomalous  component  (see  Section  2.0).  During  the 
remaining  ~  2  per  cent  of  the  time  the  particle  environment  Is  dominated  at 
low  and  moderate  energies,  by  solar  particles.  Figure  4.1  taken  from  King 
(1974),  shows  the  proton  fluence  (E  >  30  MeV)  and  the  time  of  occurence  of 
all  the  major  solar  flares  from  1S56  to  1972.  From  this  figure  we  can  see 
that  the  frequency  distribution  has  a  period  of  -  n  years  (the  solar 
cycle)  and  that  each  active  period  displays  one  anomalously  large  event. 

Apart  from  these  two  events,  the  remaining  events  seem  to  be  distributed  as 
though  the  log  (to  base  10)  of  their  sizes  was  normally  distributed.  This  Is 
called  a  log-normal  distribution  (see  Brown,  1957).  King  has  found  the 
log-normal  fluence,  F,  distributions  (Log  iof )  for  the  fluence,  F,  above 
any  energy  threshold  for  these  events.  The  means  and  standard  deviations  of 
the  distributions  for  four  thresholds  are  shown  In  Table  4.1. 

The  frequency  distribution  of  flares  is  best  described  by  the  Burrell 
distribution  (see  Burrell,  1971), 

p(n,t,N,T)  -  (n  +  N) i  ( t/T)n/[n!N!  ( 1  +  t/T)1+n+hi]  (4.1) 


where  p  Is  the  probability  that  exactly  n  flares  will  occur  during  a  time  t 
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given  that  N  flares  were  observed  curing  a  time  T.  For  the  second  active 
period  In  Figure  4.1,  we  have  N  ■  24,  T  «  7  years. 

Because  of  Its  size,  the  August  1972  event  (actually  a  series  of  four 
flares),  produced  a  fluence  at  earth  nearly  twice  that  of  all  the  other  flares 
in  this  active  perloo  combined.  King,  therefore,  treats  this  flare 
separately.  If  a  satellite  mission  is  long  enough  for  the  risk  o*  such  an 
event,  p(l,t,l,7)  In  Eq.  (4.1),  to  be  unacceptable,  then  it  must  be  Included 
in  the  particle  environment  model.  When  anomalously  large  events  are 
considered  they  always  dominate  ordinary  events  In  total  fluence  for  the 
mission. 

A  solcr  particle  event  may  last  several  hours  or  days,  during  this  time 
the  flux  varies  enormously.  Besides  knowing  the  integral  fluence  for  the 
event,  It  Is  useful  to  know  the  maximum  flux.  Using  data  provided  by  King 
(1974)  we  have  a  log-normal  distribution  to  the  peak  proton  fluxes  above  three 
energy  thresholds.  The  values  of  the  means  and  standard  deviations  are  shown 
in  Table  4.1: 


TABLE  4.1  The  Parameters  of  Log-normal  Distributions  for  Ordinary 
Solar  Flares.  Parameters  are  showr  for:  (a)  the  ^Integral  omni-directional 
fluence  for  the  entire  solar  event  in  proton s/cm^  and  (b)  the  peak  omni- 


directional 

flux  in  protons/cm*1  sec. 

E  >  10  MeV 

E  >  30  MeV 

E  >  60 

MeV 

E  >  100  MeV 

Integral  Fluence 
mean  +  o 

8.27  +  .59 

7.28  +  .75 

6.63  + 

.95 

5.77  +  1.24 

Peak  Flux 
mean  +  o 

3.27  +  .£4 

2.37  +  .82 

1.88  + 

.78 

4.2  Solar  Proton  Spectra 

From  Table  4.1  we  can  see  that  the  ordinary  solar  flare  proton  spectra  vary 
enormously  in  amplitude  and  spectral  shape.  The  mean  log  fluence  for  ordinary 
events  has  been  fit  to  give  a  differential  spectrum  of  the  form. : 

Fmean  “  3,3  x  10s(e'E/,2°*2  +  307e'E</3)  protons/cmzster.MeV  (4.2) 

for  E  >  10  MeV/u. 

This  is  the  typical  spectrum  of  particles  that  we  expect  to  arrive  In  the 
1 nterplanetary  medium  near  earth  as  the  result  of  an  ordinary  flare.  Integrated  over 
the  period  of  the  flare.  A  "worst  case"  spectrum  can  be  obtained  by  fitting  the 
mean  +  1  o  for  each  threshold  in  Table  1.  If  the  log-normal  distributions  were 
uncorrel a ted,  this  would  produce  a  case  which  would  be  exceeded  with  a  probability 
of  only  0.014.  The  distributions  are,  of  course,  correlated  to  some  degree  so  the 
probability  Is  somewhat  higher,  but  not  larger  than  0.34.  An  extensive  study  of  the 
correlations  between  the  spectra  would  be  required  to  determine  this  probability. 

For  now,  we  will  assume  that  the  spectra  have  two  Independent  parameters  and  take 
the  probability  to  be  0.12. 


38 


We  have  also  fit  the  "worst  case"  spectrum,  as  shown  below: 
Fworst  *  7.6  x  lC5(e-t'/30  +  165  e-E/4.0)  protons/cm^ster.MeV 


(4.31 


for  E  10  MeV.  As  eq.  (4.t)  above,  this  Is  the  most  Intense  spectrum  (with 
a  X)  per  cent  confidence  level)  that  we  expect  to  find  In  the  Interplanetary 
medium  near  earth  as  the  result  of  an  ordinary  flare.  Both  Fmpan  anH 
F worst  are  shown  In  Figure  4.2. 

The  same  fits  have  been  done  for  the  peak  flux  distributions  (l.e.,  the 
part de  flux  at  the  peak  of  the  flare's  intensity),  they  are  shown  below  and 
in  Figure  4.3  for  toe  mean  (typical)  and  worst  (90  per  cent)  cases. 


fmean  3  1 . 95 ( e“E/27 . 5  +  l 73e “E/4 )  protons/cm2ster.sec.MeV  (4.4) 

and 

^worst  -  17. l[e“E/24.5  +  63.6e-E/4]  protons/cm2ster .sec. MeV  (4.5) 

These  are  the  mean  and  worst-case  spectra  to  be  expected  during  the  most 
Intense  part  of  the  flare. 

For  anomalously  large  events,  King  (1974)  suggests  the  event 
time-integrated  spectrum  In  the  Interplanetary  medium  near  earth: 


Fa  =  2.37  x  10?exp[(30  -El/26.5]  protons/cm2ster .MeV  (4.C) 


This  spectrum  is  compared  with  eqs.  (4.2)  and  (4.3)  In  Figure  4.2. 

Figure  4.4  taken  from  Lockwood  et  al .  (1575),  shows  the  existing  data 
on  the  peak  flux  spectrum  for  the  August  4,  1972  event.  It  gives  a  feeling 
for  the  kind  of  uncertainty  that  exists  In  the  measurements  of  the  spectrum 
oJ  nT9r  f]are‘  For  the  purpose  of  this  model,  we  recommend  extending  the 
20:00  UT  Explorer  41  spectrum  to  150  MeV  and  matching  it  to  a  P-9  power  law 
at  that  point, 

That  Is: 

f  «  dP  x  g.3  x  iq5  e-P/0.10  protons/cm2ster.sec.  MeV  (4.7) 


for  E  <  150  MeV. 
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arid, 

f  =.  dP  x  i7,g  p-9  protons/cm2ster.sec .MeV  (4.8) 

dt 

for  £  >  ISO  MeV, 
whore 

P  “  C ( E/1000 )2  +  1.86  x  10-3e]1/2  (4.9) 

with  £  In  MeV.  This  gives  the  particle  energy  spectrum  to  be  expected  in  the 
interplanetary  medium,  near  earth,  during  the  most  Intense  part  of  an 
anomalously  large  flare.  It  is  compared  to  eqs.  (4.4)  and  (4.5)  in  Figure  4.3 


4.3  Solar  Energetic  Particle  Composition 

The  elemental  composition  of  particles  from  solar  flares  is  highly 
variable  showing,  in  some  cases,  enormous  enhancements  in  heavy  elements. 

The  data  on  several  large  solar  flares  between  October  30,  1973  and  December 
1,  1977  have  recently  been  surveyed  by  Mason  et  al.  (1980).  The  average 
composition  they  found  at  ~  1  MeV/u  for  37  days  during  major  flares  Is 
given  in  Table  4.2,  normalized  to  hydrogen. 


TABLE  4.2  Solar  Energetic  Particle  Composition* 


Element 

N 

He 

C 

0 

Me 

Mg 

SI 

S-Ca 

Cr-Ni 

Mean 

1 

.022 

1.64 

3.24 

5.16 

4.8£ 

3.85 

2.55 

4.45 

Mean  +  lo 

1 

.031 

3.54 

8.74 

1.64 

1.54 

1.34 

9.55 

2.04 

Mean  -  lo 

1 

.020 

1.4« 

2.34 

3.05 

2.85 

1.85 

l.l5 

1.4S 

Richest  day 

1 

.074 

.83 

3.23 

(j.84 

7.74 

8.14 

4.84 

l.l3 

"Notation  in  this 
.0C0Q11 

table 

has 

been  compressed 

1.15  means  1.1  x 

10-5  or 

These  results  are  consistent  with  a  more  limited  survey  carried  out  by 
Webber  (1975)  for  E  >  20  MeV. 

Mason  et  al.  find  that:  (1)  the  average  composition  does  not  depend 
strongly  on  particle  energy  or  flare  size;  (2)  all  extreme  examples  of 
composition  anomalies  are  for  small  flares. 

Comparing  Table  4.2  with  the  cosmic  ray  composition  described  in 
Section  2.0,  we  see  that  solar  flares  have  a  H/(C  +  0)  ratio  about  ten 
times  larger  than  cosmic  rays.  The  He/(C  +  0)  ratio  Is  ~  45  In  solar 
flares  compared  to  21  in  cosmic  rays  while  the  Fe/(C  +  0)  ratio  appears  to 
be  ~  1,5  times  larger  in  flares  than  cosmic  rays.  This  Is  due  to 
spallation;  the  Fe/(C  +  0)  ratio  at  cosmic  ray  sources  Is  larger  than  in 
flares.  The  elements  Li,  Be,  B  and  the  odd  2  elements  above  nitrogen  are 
much  less  abundant  In  solar  flare  particles  than  cosmic  rays.  This  is 
because  spallation  has  filled  in  these  nuclei  in  cosmic  rays. 
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Table  4.2  also  gives  a  feeling  for  the  variability  of  the  dally 
average  composition  during  enhanced  periods  associated  with  large  flares. 
The  mean  +  1  a  line  gives  the  limit  on  heavy  Ion  richness  that  Is  exceeded 
by  only  one  In  every  six  dally  averages.  Correspondingly,  the  mean  -1  a 
line  gives  composition  limit  for  heavy  Ion  poor  flares;  only  one  In  six 
days  were  poorer  In  heavy  Ions.  It  should  be  noted  that  the  error  bounds 
are  not  equidistant  on  either  side  ol  the  mean,  the  distribution  has  a  much 
larger  wing  to  the  heavy-ion-rich  side  and  the  distribution  Is  not  gausslan 
shaped.  To  give  a  feel  for  the  worst  case,  the  last  line  of  the  table 
gives  the  results  for  the  richest  day  of  the  37  days  Included  In  the 
survey.  Mason  et  al.  also  show  that  these  enrichments  are  highly 
correlated  from  element  to  element. 

The  survey  of  Mason  et  al.  covers  11  periods  of  flare  activity.  Cook 
et  al .  (1980)  report  results  for  Z  >  2  measured  In  the  4.6  to  8.7  MeV/u 
energy  range  from  four  flares  In  1978.  Their  results  show  that  one  flare, 
April  29,  1978  was  unusually  rich  in  all  heavy  elements.  The  He/(C  +  0) 
ratio  was  10,  a  factor  of  4.5  richer  than  the  average  flare  and  a  factor  of 
2  richer  than  cosmic  rays. 

The  available  data  on  solar  flare  composition  Is  still  quite  limited 
at  low  energies  and  very  sparse  at  the  higher  energies  of  Interest  here. 

The  variations  In  composition  from  flare  to  flare  are  large  and  distributed 
In  a  broad  non-gausslan  manner.  The  result  Is  that  the  uncertainty  In  the 
flux  of  any  elemental  species  Is  due  almost  as  much  to  the  variations  In 
the  composition  of  energetic  particles  as  to  the  variations  In  flare  sizes. 

To  obtain  a  worst  case  composition  for  a  given  confidence  level,  we 
have  chosen  to  treat  the  composition  as  though  It  were  normally  distributed 
with  different  standard  deviations  above  and  below  the  mean,  using  the 
standard  deviations  from  Table  4.2,  1 nterpolatlng  to  neighboring  elements 
as  needed.  The  abundances  of  the  elements  P,  Cl,  K,  T1,  Mn,  and  Co  were 
estimated  from  the  solar  system  abundances  compiled  by  Cameron  (1980). 

Table  4.3  gives  our  recommendation  for  the  mean  composition  and  the 
worst  case  composition  at  the  90  per  cent  confidence  level  (l.e.,  there  is 
only  one  chance  In  ter  of  having  a  richer  flare).  Elements  with  mean 
relative  abundances  below  10-7  *ere  treated  as  absent  In  the  composition. 
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TABLE  4.3  Mean  and  Worst  Case  Compositions 


Mean  Case 

Worst 

Case 

Mean  Case 

Worst  Case 

H 

1 

1 

p 

2.3 

X  10"7 

1.1 

x  10“ 6 

He 

2.2 

X 

10‘2 

3.3 

X 

10"2 

s 

1.8 

x  10-5 

8.4 

x  10"5 

LI 

0 

0 

Cl 

1.7 

x  10"7 

8  x 

10“7 

Be 

0 

0 

Ar 

3.9 

x  10“6 

1.8 

x  10"5 

B 

0 

0 

K 

1.3 

x  10“7 

6  x 

10“7 

C 

1.6 

X 

10“4 

4.C 

X 

10“4 

Ca 

2.3 

x  10“6 

1  X 

10“5 

N 

2.8 

X 

10'5 

1.1 

X 

10"4 

Sc 

0 

0 

0 

3.2 

X 

10"4 

1.0 

X 

10’3 

T1 

1  X 

10"? 

5  x 

10“7 

F 

0 

0 

V 

0 

0 

tie 

5.1 

X 

10-5 

1.9 

X 

10'4 

Cr 

5.7 

X 

o 

1 

3.2 

x  10“e 

Na 

1.6 

X 

10'6 

6.1 

X 

10"6 

Mn 

4.2 

x  10“7 

2.3 

x  10'6 

Mg 

4.8 

X 

10-5 

1.8 

X 

10“4 

Fe 

4.1 

x  10“£ 

2.3 

x  10“4 

Al 

3.6 

X 

10'6 

1.4 

X 

10-5 

Co 

1  X 

10’7 

5.5 

X 

o 

1 

-V4 

SI 

3.8 

X 

icr5 

1.6 

X 

10"4 

N1 

2.2 

to 

o 

r  -t 

X 

1.2 

x  10“5 

Z>28 

0 

0 

Besides  the  atomic  nuclei,  flares  accelerate  electrons.  As  discussed  In 
Section  2.4,  electrons  will  only  be  Important  If  they  are  overwhelmingly 
more  abundant  than  nuclei  at  energies  greater  than  20  MeV.  Samaty  et  al. 
(1980)  argue,  based  on  the  existing  data,  that  the  electron  to  proton  ratio 
at  energies  greater  than  10  MeV  Is  ~  10”2,  clearly  solar  flare  electrons 
are  not  a  problem  and  will  be  Ignored  In  this  model. 

4.4  Recommendations 


Large  solar  flares  are  transient  phenomena,  contributing  to  the  particle 
environment  only  about  2  per  cent  of  the  time.  If  a  few  per  cent  operating 
time  can  be  lost  during  flares,  then  their  contribution  to  the  particle 
environment  can  be  Ignored  In  formulating  a  worst  case  to  be  withstood. 

If  flares  must  be  considered,  then  we  recommend  the  following  procedure: 
(1)  determine  whether  there  Is  an  unacceptable  risk  of  an  anomalously  large 
event  from  Eq.  4.1  with  N  *  1,  T  ■  7;  (2)  If  anomalously  large  events  must 
be  considered,  assume  the  event  Integral  spectrum,  Eq.  4.6,  and  peak  flux 
spectrum  Eqs.  4.7,  4.8,  and  4.9  for  protons;  (3)  If  anomalously  large  events 
are  unlikely,  use  Eq.  4.1  to  estimate  the  number  of  ordinary  events  to  be 
expected;  (4)  Eqs.  4.3  and  4.5  give  the  worst  case  event  fluence  and  peak 
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flux  proton  spectra  to  be  considered;  (6)  In  any  case,  Table  4.3  should  be 
used  to  obtain  the  worst  case  composition.  Multiplying  the  proton  spectra 
by  the  numbers  In  this  table  will  give  estimates  of  the  spectra  for  the 
other  elements. 

When  spacecraft  Inside  the  magnetosphere  are  studied,  the  geomagnetic 
cutoff  should  be  taken  Into  account  (see  Section  5.0),  using  the  model  for  a 
ulsturbed  magnetosphere. 

Some  thought  should  be  given  to  the  problem  of  modulating  the  peak  flux 
spectrum  with  the  geomagnetic  cutoff  transmittance  function.  If  the  orbital 
period  Is  ~  100  minutes  then  the  flare  peak  will  be  (at  least  partially) 
averaged  over  the  orbit.  For  long-period  orbits,  the  flare  peak  will  not  be 
averaged  and  a  position  for  the  satellite  must  be  assumed. 
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ENERGY  (M*v) 

Fig,  4.3  —  The  peak  proton  flux  differential  energy  ipectra  for:  F_ ,  a  typical 
ordinary  event;  Fw,  a  worst-caw  ordinary  event  (80  percent  confluence  level); 
and  F,  anomalously  large  event. 
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Fi^.  4.4  —  Differential  lolnr  proton  magnetic  rifeidity  spectra 
during  the  August  4  event  (from  Lockwood  ot  ul.  1975). 


The  earth's  magnetic  field  must  be  penetrated  by  cosmic  rays  in  order 
for  them  to  reach  a  spacecraft  in  earth  orbit.  The  number  of  magnetic  field 
lines  a  cosmic  ray  must  cross  to  reach  a  given  point  within  the 
magnetosphere  approximately  determines  the  minimum  energy  It  must  possess. 

To  cross  more  magnetic  Held  lines  more  energy  will  be  required.  This 
penetrating  ability  Is  determined  uniquely  by  the  cosmic  ray's  momentum 
divided  by  Its  charge.  This  quantity  is  called  magnetic  rigidity  (see 
Rossi,  1964,  Appendix  F).  To  penetrate  the  earth's  magnetic  field,  a 
particle  must  have  sufficient  magnetic  rigidity  (momentum  per  unit  charge) 
to  avoid  being  turned  away.  There  is  a  minimum  magnetic  rigidity  a  cosmic 
ray  must  possess  to  arrive  from  a  given  direction  at  a  given  point  in  the 
magnetosphere.  Regions  In  the  outer  mignetosphere  and  near  the  poles  can  be 
reached  at  much  lower  magnetic  rigidities  than  are  required  to  reach  points 
near  the  earth's  equator.  In  general,  for  each  point  in  the  magnetosphere 
and  for  ,ach  direction  from  that  point,  there  exists  a  magnetic  rigidity 
below  wUch  cosmic  rays  cannot  arrive.  This  value  Is  the  geomagnetic 
cutoff.  For  magnetic  rigidities  above  this  value,  cosmic  rays  arrive 
freely,  as  though  no  magnetic  field  were  present. 


5.1  Methods  for  Computing  the  Cutoff 


The  geomagnetic  cutoff  was  first  calculated  by  C.  Stomter  (1930),  using 
a  dipole  approximation  for  the  earth's  magnetic  field.  He  showed  that  the 
cutoff  rigidity  at  the  earth's  surface  Is  given  by: 


P  *  —■  Ll-d-cosy  cos3aW/[co.y  cos  x]2 
r* 

for  positively  charged  particles,  where 


(5.1) 


P  ■  magnetic  rigidity  In  CeV/ec, 

r  -  radial  distance  from  the  dipole  center  In  earth  radii 
\  »  latitude  In  dipole  coordinates  and, 
y  *  the  angle  which  trajectory  makes  with  magnetic  west. 


The  magnetic  rigidity,  P,  Is  related  to  the  particles  energy  by: 


E  -  (wjj  +  P2Z2/A2)1/2  -  M0  (5.2) 

where  £  Is  the  kinetic  energy  in  GeV/u  (1  GeV/u  *  1000  MeV/u), 

P  Is  the  magnetic  rigidity  In  GeV/ec 
A  Is  the  particle's  mass  in  atnu 
2  Is  the  particle's  charge  and 

M0  -  0.931  GeV 
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Stormer' s  theory  does  not  account  for  the  presence  of  the  solid  earth, 
so  In  some  directions  at  each  location,  this  theory  predicts  a  cutoff  that 
Is  lower  than  the  actual  cutoff.  The  problem  of  the  earth's  shadow  was 
first  addressed  by  Vallarta  (1948),  again  In  the  context  of  the  dipole 
model.  Vallarta  showed  that  there  existed  a  range  of  magnetic  rigidities 
above  the  Stormer  cutoff  where  the  earth's  shadow  casts  a  broken  pattern  of 
allowed  and  forbidden  bands  of  magnetic  rigidity.  There  Is  the  penumbra! 
shadow  of  the  earth.  The  width  of  the  shadow  varies  from  10  per  cent  to  100 
per  cent  above  the  Stormer  cutoff  at  the  earth's  surface  for  zenith  angles 
<  45",  For  larger  zenith  angles,  the  effect  Increases  as  the  arrival 
cflrectlon  approaches  the  horizon.  The  density  of  the  penumbral  shadow  Is 
also  highly  variable. 

While  these  early  Investigations  revealed  the  basic  features  of  the 
geomagnetic  cutoff,  they  were  limited  in  their  accuracy  because  they 
depended  upon  a  dipole  field  model  fit  to  the  true  geomagnetic  field. 

Shea  and  Smart  (1975)  have  made  detailed  calculations  of  the  geomagnetic 
cutoff  using  a  realistic  field  model  (IGRF-196E,  1969).  This  model 
describes  the  earth's  field  In  an  80  term  spherical  harmonic  expansion.  For 
this  reason,  no  analytic  solution  exists  to  the  Stormer  problem.  Shea  and 
Smart  have  calculated  the  geomagnetic  cutoff  by  detailed  ray  tracing 
backwards  from  the  point  of  Interest,  to  determine  If  the  trajectory  leads 
back  to  outer  space.  Figure  6.1  (taken  from  Lund,  1980)  shows  such  a 
computed  trajectory  for  a  cosmic  ray  whose  magnetic  rigidity  Is  barely 
adequate  to  bring  It  Into  the  atmosphere  over  Slnklang  Province  In  Western 
China.  Such  a  calculation  must  be  carried  out  for  each  point,  each 
direction  and  each  magnetic  rigidity  of  interest  using  a  high  speed 
computer.  This  technique  has  provided  very  realistic  maps  of  the 
yeoniagnetlc  cutoff  as  well  as  detailed  examinations  of  the  penumbral 
shadow.  The  authors  define,  tor  each  point  ano  direction:  (1)  a  cutoff, 
below  which  no  open  trajectories  are  found,  the  Stormer  cutoff;  (2)  a  cutoff 
above  which  there  are  no  closed  trajectories,  the  main  cutoff  and;  (3)  the 
effective  cutoff,  a  value  between  the  main  and  Stormer  cutoffs  weighted 
according  to  the  density  of  the  penumbral  shadow. 

The  size  of  the  computational  task  requirea  to  employ  this  technique  Is 
such  that  It  can  only  be  used  for  selected  sites  and  directions.  Shea  and 
Smart  11975)  provide  calculations  of  the  vertical  geomagnetic  cutoff  at  an 
altitude  of  10  Km  on  a  world-wide  grid  of  points  spaced  apart  15*  In 
longitude  by  5*  In  latitude. 

To  obtain  the  cutoff  In  other  directions  and  at  other  points  on  earth, 
Shea  et  al  (1973)  recommend  that  the  Stormer  theory  be  used  to  Interpolate 
the  computed  vertical  cutoffs  from  the  world  grid.  These  authors  show 
excellent  agreement  between  transmittance  functions  for  cosmic  ray 
experiments  at  Palestine,  Texas,  calculated  both  by  ray  tracing  and  by 
Stormer  theory  Interpolation, 

Smart  and  Shea  (1977)  show  that  using  eq.  5.1  with  the  "best  fit" 
eccentric  dipole  model,  which  they  suggest,  It  is  possible  to  Interpolate 
from  the  world  grid  (Shea  and  Smart,  1976)  In  three  dimensions  provided  that 
radial  Interpol atlons  are  to  altitudes  small  compared  to  an  earth  radius. 


Heinrich  and  Spill  (1979)  have  used  these  Ideas  to  calculate  the 
vertlca1  geomagnetic  cutoffs  at  various  altitudes.  The  authors  also 
calculated  the  vertical  transmittance  functions  for  several  223  Km  circular 
satellite  orbits,  see  Figure  5.2.  The  technique  of  Heinrich  and  Spill 
(1979)  appears  to  be  the  best  for  calculating  the  transmittance  function  for 
an  orbit. 

The  transmittance  function  calculated  by  Heinrich  and  Spill  needs  to  be 
extendeo  to  Include  cutoffs  from  all  directions,  Including  vertical.  This 
can  be  done  [as  Smart  and  Shea  (1977)  have  suggested]  with  Stormer  theory, 
(see  eq.  3.1)  extrapolating  from  the  400  Km  altitude  world  grid  recently 
prepared  by  these  authors  (Smart,  1980).  To  estimate  the  cutoff  at  higher 
altitudes,  It  Is  necessary  to  scale  by  1/r2  from  these  calculations  since 
none  exist  at  higher  altitudes. 

There  remains  one  problem,  how  to  account  for  the  earth's  umbra! 
shadow.  On  the  earth's  surface  this  Is  simple,  cosmic  rays  can  arrive  from 
above,  not  below.  At  satellite  altitudes  the  problem  Is  not  so  simple,  for 
the  highest  energies,  the  portion  of  the  geometry  factor  that  Is  occulted 
falls  off  with  altitude,  h,  as 

Omega  •  2n  {  1  -  L(Rfi  +  hj2  -  Rfi2]  1/2/(Sfe  +  h)  \  (6.3) 


where  Rc  Is  the  earth's  radius.  At  lower  magnetic  rigidities,  the  earth's 
umbral  shadow  Is  distorted  by  the  earth's  field  and  swept  off  to  an  easterly 
direction  so  that  particles  may  arrive  below  the  optical  horizon  In  the 
west.  This  distortion  Increases  at  lower  rigidities  as  the  cutoff  Is 
approached.  Besides  the  change  of  direction  of  earth  occultatlon  at  low 
rigidities,  the  occulting  solid  angle  also  falls  off  more  rapidly  with 
altitude  than  describee  by  eq.  5.3  (Smart,  1980).  The  details  of  how  the 
earth's  umbral  shadow  changes  with  altitude  and  rigidity  are  unknown;  Smart 
(1980)  has  suggested  that  the  problem  might  be  solved  by  ray  tracing  at  a 
range  of  altitudes  and  rigidities. 

S.2  The  Effect  of  N agnetlc  Storms 

So  far  In  this  dlscusslun  we  have  only  dealt  with  the  quiescent 
magnetosphere.  When  a  solar  flare  occurs,  It  usually  causes  a  magnetic 
storm  at  aarth.  This  storm  disrupts  the  magnetosphere  altering  the 
geomagnetic  cutoff,  usually  depressing  It.  Figure  5.3  shows  the  fractional 
depression  aP/P0  In  geomagnetic  cutoff  as  a  function  of  quiescent  cutoff 
pQ,  In  two  magnetic  storms;  Nov.  15,  1950  (Webber,  1962)  and  April  1,  1973 
(Debrunner  and  Flucklaer,  1977).  The  effect  seems  to  be  the  result  of  ring 
currents  Induced  by  the  sudden  commencement  of  the  storm  (see  Flucklger  et 
al .  1979,  Dorman,  1974).  These  currents  reduce  the  equatorial  magnetic 
field  by  -  lO^y  (where  y  •*  lO"5  gauss),  allowing  penetration  to  any 
given  point  In  the  magnetosphere  by  lower  energy  cosmic  rays  than  Is 
normally  possible. 
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As  a  model  for  this  effect,  we  recommend  the  function  shown  an  the  solid 
curve  In  Figure  5.3.  Specifically, 


Ap/p0  -  .54  exp(-P0/2.S  GV)  with  Pst0rm  “  P0  "  Ap 


(5.4) 


5.3  Recommended  Procedure 


Based  on  the  foregoing  discussion,  we  recommend  for: 

1)  The  quiescent  cutoff.  Use  of  Shea  and  Smart's  400  Km  world  yrld  of 
effective  cutoffs,  Interpolating  to  other  latitudes,  longitudes  and 
altitudes  using  Eq.  1  with  the  eccentric  dipole  recommended  by  Smart  and 
Shea  ( 1377) .  The  earth's  shadow  should  be  taken  as  shrinking  with  altitude 
according  to  Eq.  2.  The  transmittance  function  Is  best  computed  by  the 
orbit  Integration  technique  of  Heinrich  and  Spill  (1979). 

2)  During  solar  flares.  Calculate  the  quiescent  cutoff  as  described 
above  and  use  Eq.  5.4  to  obtain  the  depressed  cutoff  for  flare  conditions. 
There  are  several  simple  but  less  accurate  alternatives  to  the  above 
procedure.  The  simplest  is  to  use  the  rule  of  thumb  (Smart  and  Shea,  1967). 

P0  *  15.96/L2-006 

where  L  Is  the  Mcllwaln  L  parameter  (Mcllwaln,  1961)  and  P0  is  the 
vertical  cutoff  In  GV.  P0  may  then  be  taken  as  the  average  cutoff  for  the 
directions  not  occulted  by  the  earth  (see  Eq.  5.3).  A  better  technique  Is 
to  use  the  1/L2  rule  to  extrapolate  the  400  Km  world  grla.  This 
essentially  Is  what  was  done  by  Heinrich  and  Spill  (1979). 

The  principal  sources  of  error  are  the  uncertainty  in  the  cutoff  and  the 
darkness  of  the  penumbral  shadow  at  points  reached  by  extrapolation.  This 
probably  leads  to  no  more  than  a  +  20  per  cent  error  in  the  cutoff.  A 

second  source  of  error  Is  the  size  we  assume  for  the  umbral  shadow;  this 

could  be  overestimated  by  a  factor  of  2.  This  error  would  be  about  the  same 

as  a  +  10  per  cent  error  In  magnetic  rigidity  at  400  Km  and  grow  smaller  at 

higher  altitudes.  The  actual  error  In  particle  flux  resulting  from  these 
errors  In  rigidity  depends  on  where  the  cutoff  comes  on  the  differential 
energy  spectrum,  but  even  at  high  cutoff  values  these  errors  translate  Into 
no  more  than  a  +  30  per  cent  error  In  the  particle  flux. 

The  uncertainty*  In  the  actual  value  of  the  geomagnetic  cutoff  may  be 
»  50  per  cent  at  low  cutoffs  during  a  magnetic  storm.  This  translates 
Into  only  a  +  50  per  cent  error  In  the  flux,  and  since  this  condition  Is 
transient,  IT  Is  only  Important  for  solar  flares.  The  flux  uncertainty 
during  a  solar  flare  is  >>  50  per  cent. 
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Fig.  5,2  —  The  geomagnetic  transmittance  function  for  22C  Km  circular  orbits 
of  various  inclinations.  The  symbol  R  is  used  here  for  magnetic  rigidity  taken 
from  Heinrici;  and  Spill  (1979). 
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P„  IN  GV 

Fig.  5.3  —  The  fractional  depression  of  the  geomagnetic  cutoff,  AP/P0,  as 
a  function  of  quiescent  cutoff,  P0.  The  data  are  for  magnetic  storms,  Nov. 
15,  1960,  and  April  1,  1973.  The  solid  curve  is  the  fractional  cutoff  de¬ 
pression  function  we  recommend  using  to  describe  the  geomagnetic  cutoff 
for  major  solar  flares. 
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6.0  Particles  from  the  Magnetosphere 
6.1  Protons 


Protons  are  the  most  abundant  particles  in  the  magnetosphere  that  can 
easily  produce  soft  upsets.  The  trapped  preton  environment  has  been 
comprehensively  studied  and  is  well  descrioed  by  the  computer  model ,AP-8 
(Sawyer  and  Vette,  1976).  We  recommend  that  this  model  be  used  for  the 
trapped  proton  environment.  In  addition  to  protons,  the  magnetosphere 
contains  helium  nuclei  as  well  as  heavier  nuclei,  especially  carbon, 
nitrogen  and  oxygen.  While  these  nuclei  are  less  abundant  than  protons 
they  are  much  more  effective  i n  producing  soft  upsets. 

6.2  Alpha  Particles 

Helium  nuclei  (mostly  alpha  particles)  have  been  detected  throughout 
the  magnetosphere.  The  principal  source  of  these  nuclei  appears  to  be  the 
solar  wind  (Blake,  1973  and  Hovestadt  et  al.  1978).  The  solar  wind 
particles  are  transported  down  into  the  magnetosphere  and  accelerated  by 
radial  oiffuslon.  This  process  was  described  theoretically  by  Cornwall 
(1572)  and  has  recently  been  shown  to  describe  well  the  helium  Ion 
population  In  the  magnetosphere  (Spjeldnlk  and  Fritz,  1978,  and  Fritz  and 
Spjeldnik,  1979).  The  bulk  of  the  helium  nuclei  arc,  however,  at  energies 
too  low  to  penetrate  the  walls  of  the  spacecraft. 

H.e  only  measurements  of  geomagnetlcally  trapped  alpha  particles  at 
energies  above  ~  2  MeV/u  have  been  reported  by  Rubin  et  al  (1977)  and 
Panasyuk  et  al.  (1977).  The  results  of  Rubin  et  al  on  the  helium  to 
proton  ratio  are  shown  In  Figure  6.1.  Also  shown  for  comparison  are  the 
low  energy  measurements  of  Fennell  et  al .  (1974)  and  Blake  et  al  (1973). 
The  data  of  Rubin  et  al  cover  the  range  of  1.85  <  L  <  2.65  at  low 
altitudes,  l.e.  3  <_  B/Bq  _<  8.5.  The  data  point  of  Fennell  et  al  is  at 
2.6  <  L  2.7,  1.3  _<  B/Bo  <  2.3,  while  that  of  Blake  et  al.  Is  on  the 
geomagnetic  equator  at  L  ■  1.95.  The  results  of  Panasyuk  et  al.  (1977) 
cove-'  the  range  from  2.0  <  L  <  2.55  and  1.4  _<  B/B^  <  2.35.  Their  data 
cover  the  broad  energy  band  4.25  <  E  <  15  MeV/u.  For  this  range,  they 
report  5  x  10“^  _<  «/p  <  1.5  x  10"^.  These  results  are  consistent  with 
those  shown  in  Figure  6.1. 


These  results  may  be  taken  as  typical  of  quiet  periods  In  the  Inner 
zone,  although  there  was  one  magnetic  storm  during  the  period  (Rubin  et 
al.  1977). 

The  declining  trend  in  the  «/p  ratio  In  Figure  6.1  can  be  expected  to 
continue  co  higher  energies.  This  Is  because  the  protons  are  principally 
from  cosmic  ray  albedo  neutron  decay  (CRAND)  while  the  helium  nuclei  have 
diffused  In  from  some  external  source,  presumably,  the  solar  wind. 
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Rubin  et  al  (1977)  show  that  the  a/p  ratio  varies  about  a  factor  of  2 
between  l  values  of  1.85  and  2.65  with  a  peak  at  L  ■  2.45.  For 
purposes  of  this  particle  environment  model,  we  have  chosen  to  Ignore  the 
L  value  dependence  of  the  a/p  ratio  and  adopt  the  dashed  curve  In  Figure 
5.1  as  the  energy  dependent  a/p  ratio  for  L  <  2.6.  Specifically, 

a/p  -  2.5  x  10"3;  £  <  5  MeV/u  (6.1) 

a/p  -  8.3  x  10"3  exp( -E/4.15);  5  MeV/u  <  E  _<  10  MeV/u  (6.2) 

and 

a/p  *  0.957  exp ( -E/1 . 4 ) ;  E  >  10  f-ieV/u  (6.3) 


There  are  no  quiet  time  measurements  of  helium  nuclei  above  10  MeV/u  at 
L  >  2.65  so  we  must  use  the  data  at  lower  energies.  Fritz  and  Spjeldnlk 
(1579)  have  computed  a/p  ratios  at  the  same  energy  per  amu  that  best  fit  the 
data  over  a  range  of  L  values.  Based  on  their  results,  we  adopt  the 
conservative  value  of  a/p  ■  2.5  x  10‘3  for  L  >  2.5.  While  the 
calculations  of  Fritz  and  Spjeldnlk  extend  only  to  2  MeV/u,  we  have  made  the 
assumption  that  this  ratio  Is  the  same  at  energies  of  2  MeV/u  or  higher. 

It  should  be  noted  that  above  L  ■  2.5,  the  proton  flux  with  E  >  5  MeV 
falls  off  very  rapidly  so,  In  practice,  most  of  the  helium  nuclei  with  E  > 

15  MeV/u  will  still  be  found  below  L  ■  2.5. 

There  have  been  a  few  reports  of  a/p  ratio  measurements  during  solar 
flares  and  magnetic  storms.  These  results  differ  considerably  from  the 
quiet  time  model  presented  above.  Yerzarlu  (1973)  reported  results 
following  the  solar  flare  of  March  C,  1970.  This  flare  had  a  proton  fluence 
(>  30  MeVj  of  1.3  x  106/cm‘,  (compared  to  8.1  x  109/cm2  for  the 
August  1972  flares).  Yerzarlu  reports  a  proton  flux  increase  of  -  10,  for 
a  total  «  flux  Increase  of  -  50.  This  condition  decayed  back  to  pre-flare 
levels  over  a  few  days. 

Transient  events  of  this  type  usually  disturb  the  outer  magnetosphere 
much  more  than  the  Inner  zone.  Scboler  et  al.  (1979)  describe 
magnetospherlc  conditions  during  a  large  geomagnetic  storm  (Dst  «  -  230) 
associated  with  a  solar  flare.  They  report  that  the  a/p  ratio  Increased  by 
a  factor  of  ~  7  at  L  ■  2.48  during  the  storm  main  phase.  The  proton  flux 
appears  to  have  Increased  by  -  10  at  the  same  time,  leading  to  an  a  flux 
Increase  -  70.  The  authors  also  report  precipitating  particles,  including 
heavy  Ions,  down  to  L  -  2.7. 

The  enhanced  a/p  ratios  reported  by  Yerzarlu  and  Scholer  et  al .  are  2.3 
x  10-3  and  2.5  x  10“3  respectively.  These  seem  to  be  covered  by  our 
conservative  choice  of  a/p  ratio  for  this  model. 

In  addition  to  the  short  term  (few  day)  enhancements  cited  above,  Van 
Allen  and  Randall  (1571)  report  evidence  that  solar  flare  u  particles  may  be 
durably  trapped  In  the  magnetosphere.  They  report  on  enhancement  of  -  40 
In  the  0.2  to  2  MeV/u  a  flux  at  3.0  <  L  _<  3.5.  The  enhancement  decayed  away 
with  a  time  constant  of  -  45  days.  The  proton  flux  at  the  same  time 
increased  by  ~  4,  implying  an  a/p  ratio  Increase  of  “  10.  The  authors 
also  examined  other  a  flux  Increases  in  the  Innerplanetary  medium  and  found 
one  additional  example  cf  durable  trapping  but  several  cases  when  a's  were 
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not  durably  trapped.  They  concluded  that  special  conditions  were  required 
for  a's  to  become  durably  trapped  in  the  magnetosphere.  Spjeldnik  and  Fritz 
(1981a)  have  recently  reported  an  additional  example  of  stable  trapping 
where  the  particles  were  presumably  from  solar  flares.  These  authors  report 
an  a  flux  Increase  of  -  30  at  L  ■  2.5  for  E  ~  0.5  toeV/u  and  a  decay  time 
>  150  days  following  the  solar  flares  of  August,  1972.  At  the  same  time,  no 
appreciable  Increase  was  seen  in  the  proton  flux. 

These  results  on  durable  trapping  of  solar  flare  particles  must  be 
extrapolated  to  the  higher  energies  of  Interest  here.  There  is  no  way  of 
knowing  how  many  a-partlcles  were  Injected  by  these  flares  at  higher 
energies.  We  can  only  guess  that  the  flux  Increases  may  have  been 
comparable.  The  decay  times  observed  at  these  low  energies  will  surely  be 
much  longer  at  energies  above  10  MeV/u  and  based  on  the  analysis  of 
Spjeldnik  and  Fritz  (1981a),  we  conclude  that  the  decay  time  at  L  »  2.5  for 
equatorially  mirroring  a's  may  be  many  years!  Decay  times  will  be  more 
rapid  for  a's  mirroring  off  the  equator  because  of  the  increased  energy  loss 
In  the  residual  atmosphere.  From  this  line  of  reasoning,  it  appears 
possible  that  the  energetic  a  flux  above  10  MeV/u  and  inside  L  *2.5  may 
originate  in  solar  flares. 

Based  on  the  scanty  data  available  at  energies  above  10  KeV/u,  it  is 
difficult  to  draw  any  conclusions.  Nevertheless,  taking  an  optimistic  view, 
as  we  said  we  would  do  in  the  Introduction,  we  will  assume  that  the  results 
of  Rubin  et  al .  (1977)  and  Panasyuk  et  al .  (1977)  represent  a  typical 
sampling  of  the  a  flux.  We  therefore  conclude  that  the  a/p  ratio  model  we 
suggest  (Eqs.  6.1,  C.2,  and  6.3)  is  sufficiently  conservative  to  describe 
the  conditions  in  the  inner  magnetosphere  most  of  the  time.  We  feel  the  a 
flux  is  best  described  by  multiplying  our  a/p  ratio  by  the  proton  flux 
predicted  by  Ap-8  (Sawyer  and  Vette,  1976). 

6.3  Heavy  Nuclei 

C,  N,  and  0  have  been  observed  in  several  experiments.  It  is  by  rto 
means  clear  that  the  particles  in  all  these  observations  were  trapped  In  the 
magnetosphere,  but  in  each  instance  the  particles  were  forbidden  direct 
access  by  the  geomagnetic  cutoff,  so  they  did  not  come  in  directly  from 
outside. 

Van  Allen  et  al .  (1970)  report  the  detection  of  heavy  nuclei,  presumably 
C,  N,  and  0  in  the  range  3.0  <  L  <  3.5,  and  0.15  <  B  <  0.2.  They  found  that 
a  ratio  CNO/u  *  2.8  x  10’3  above  ff.3  MeV/u.  Hovestaift  et  al .  (1978) 
report  substantial  fluxes  of  C,  0  and  heavier  ions  between  L  *  2.5  and  4. 
They  report  a  CNO/a  ratio  of  3  x  10’z  for  0.4  MeV/u  <  E  <  1.5  MeV/u.  The 
CNO  flux  is  dominated  by  carbon  with  C:N:Q  proportions  oT  2.7:0.28:1.  In 
contrast  with  these  results,  Blake  et  al.(1980)  report  a  CNO/a  ratio  of  6.8 
x  10'5  at  L  *  3.25,  0.35  <  B  <  0.25  and  E  >  0.25  MeV.  This  is  in  clear 
disagreement  with  the  results  of  Van  Allen  et  al .  (1970)  and  it  Is  difficult 
to  reconcile  with  the  results  of  Hovestadt  et  al.  (1978)  which  were  obtained 
near  the  geomagnetic  equator.  It  appears  that  more  measurements  will  be 
required  to  resolve  these  differences. 

As  with  the  helium  data,  these  measurements  are  at  energies  too  low  to 
cause  soft  errors.  There  is  one  measurement  at  higher  energies,  reported  by 
Mogro-Campero  (19721.  The  author  v*as  not  able  to  measure  the  CNO/He  ratio, 
but  reports  2  x  1C*5  CNO  nuclei/cm*  ster.sec.  MeV/u  at  L  ■  4  with  13  <  E 
<_  33  MsV/u  which  exceeded  the  interplanetary  flux  by  -  100.  If  these- 
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particles  are  trapped  In  the  magnetosphere  as  the  authors  argue,  they 
represent  a  large  flux.  Comparing  Mogro-Campero 1 s  CNO  flux  with  the  scanty 
data  available  on  helium  nuclei,  we  conclude  that  the  CNO/a  ratio  probably 
exceeds  1  and  that  the  CNO  flux  may  be  comparable  to  the  proton  flux  at  L  ■ 
4,  13  <_  E  <  33  MeV/u. 

In  addition  to  those  heavy  nuclei  believed  trapped,  In  the 
magnetosphere,  another  population  of  unknown  origin  has  been  observed  by 
Chan  and  Price  (1976)  arid  Biswas  et  al.  (1975).  The  results  are  based  on  a 
single  experiment  performed  outside  Skylab  from  Nov.  1973  to  Feb  1974  with  a 
stack  of  plastic  track  detectors.  For  this  reason,  the  experiment  provided 
only  the  Integral  fluence  for  the  420  Km,  50*  Skylab  orbit.  The 
orbit-averaged  flux  Is  shown  In  Figure  6.2  (Biswas  and  Durgaprasad,  1980  and 
Chan  1976).  The  figure  compares  the  measured  oxygen  spectrum  with  the 
galactic  cosmic  ray  spectrum  (GCR)  modulated  by  the  geomagnetic  cutoff  of 
the  Skylab  orbit.  Also  shown  Is  the  anomalous  component  oxygen  spectrum 
assuming  the  oxygen  Is  singly  Ionized.  As  can  be  seen,  the  flux  exceeds 
that  expectod  from  both  these  sources. 

It  Is  possible  that  the  particles  were  trapped  in  the  earth's  magnetic 
field.  If  so,  they  would  have  been  collected  only  when  the  spacecraft 
passed  through  the  South  Atlantic  anomaly  (1.3  <  L  <  1.7).  About  1.4  per 
cent  of  the  orbit  time  was  spent  In  the  anomaly.  This  implies  that  the 
trapped  particle  spectrum  would  be  the  same  as  that  shown  In  Figure  6.2  with 
the  flux  multiplied  by  -  72.  This  leads  to  a  flux  of  “  3  x  10"4  CNO 
nuclel/cmZ  ster. sec. MeV/u ,  for  13  <  E  <_  33  MeV/u  and  1.3  <  L  _<  1.7,  more 
than  an  order  of  magnitude  higher  fhan^that  reported  by  Mogro-Campero,  but 
much  less  than  the  proton  flux  at  these  L  values.  Comparing  the  orbit 
averaged  proton  flux  with  the  measured  average  oxygen  flux,  we  have 
0/p  -  1  x  10"6  at  20  MeV/u. 

A  large  number  of  explanations  have  been  suggested  for  these  results 
[see  Biswas  and  Durgaprasad  (1980)  and  Price  (1979)].  In  our  Judgment,  the 
most  likely  of  these  has  been  put  forward  by  Blake  and  Frlesen  (1977). 

These  authors  suggest  that  the  particles  of  the  anomalous  component  are 
singly  Ionized  and  that,  because  of  their  large  rigidity  at  low  energies, 
these  Ions  penetrate  deep  Into  the  magnetosphere  at  low  energies.  Some  of 
these  Ions  arrive  close  to  the  atmosphere  near  their  geomagnetic  cutoff  and 
consequently  travel  in  the  local  mirror  plane  for  trapped  particles. 

Because  of  the  large  cross  section  for  stripping  (E  >  10  MeV/u),  these 
particles  become  stripped  rapidly.  Once  stripped,  they  have  a  much  lower 
magnetic  rigidity  and  because  they  are  moving  In  the  local  mirror  plane, 
become  more  or  less  stably  trapped.  This  leads  to  a  special  trapped 
population  of  oxygen,  nitrogen,  neon  and  a  few  other  elements  which  comprise 
the  anomalous  component.  This  population  1$  also  unusual  because  Its 
equatorial  pitch  angle  distribution  would  be  double  peaked  near  36s  and 
145*,  so  that  the  flux  Is  nearly  Independent  of  B/B0|  i.e.  all  particles 
mirror  at  low  altitudes. 

The  Blake  and  Frlesen  theory  Is  further  supported  by  the  composition 
results  shown  In  Table  6.1  (Biswas  and  Durgaprasad,  1980). 
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TABLE 

6,1:  Elemental 

Composition  of  Heavy  Ions 

In  the  Magnetosphere 

Element 

Anomalous  Comp. 

Relative  Composition 

Magnetosphere  Galactic  CR 

C 

.23  +  .09 

.21  +  .019 

1.13  +  .03 

N 

.22  +  .09 

.21  +  .041 

.27  +  .02 

0 

1.0 

1.0 

1.0 

Ne 

.07  +  .04 

.08  +  .02 

.18  +  .01 

Mg 

.002  +  .002 

.006  +  .004 

.20  +  .01 

SI 

<  .02 

.004  +  .002 

.14  +  .006 

S 

— 

<  .004 

.035  +  .003 

Ar 

— 

<  .003 

,013  +  .002 

Fe  group 

— 

.05  +  .02 

.084  +  .001 

As  Table  6.1  shows,  the  heavy  Ions  In  the  magnetosphere  have  a  composition 
which  matches  the  anomalous  component  much  more  closely  than  the  galactic 
cosmic  rays,  or  for  that  matter,  the  solar  composition. 

In  addition  to  the  elements  detected  by  the  plastic  track  detectors 
outside  Skylab,  a  second  experiment  was  performed  outside  Skylab  using  a 
glass  detector.  Kretschmer  (1975)  reports  the  spectrum  of  low  energy  Iron 
group  nuclei  measured  In  this  detector.  The  spectrum  was  measured  only  up 
to  10  MeV/u.  Assuming  these  particles  are  trapped,  the  trapped  spectrum 
for  L  <  1.4  would  be: 

F  3  1.3  x  10"3  E"1'68  iron  nuclei/cm2  ster. sec. MeV/u  (6.4) 

for  1  <  E  <  10  MeV/u,  and  probably  declining  above  that  energy  more  like 
the  oxygen  spectrum,  eqs.  6.5  and  C.6.  Spjelanlk  and  Fritz  (1961)  have 
reported  Z  >  9  particles  stably  trapped  at  2  _<  L  <_  3.5  and  B/B0  <  1.6  at 
energies  ~  1  Mev/u. 

While  the  Blake  and  Frlesen  theory  explains  the  Skylab  experiment,  It 
cannot  explain  Mogro-Campero' s  QGO-5  data.  The  particles  In  this 
experiment  have  equatorial  pitch  angles  >  24*  (70  per  cent  of  theirr>  45*) 
while  the  Blake  and  Frlesen  theory  would  predict  pitch  angles  ~  16*. 

It  would  appear,  based  on  these  two  experiments,  that  we  have  a 
relatively  large  component  of  energetic  trapped  heavy  ions  In  the 
magnetosphere.  During  quiet  times  In  the  Inner  zone,  the  proton  flux  might 
be  the  dominant  cause  of  soft  errors,  but  beyond  L  -  2.5  the  heavy  lun 
flux  becomes  Increasingly  Important. 

Besides  these  "quiet  time"  measurements,  Spjeldnlk  and  Fritz  (1981b) 
have  reported  two  events  In  which  heavy  Ions  were  Injected  deep 
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into  the  magnetosphere  by  solar  flares,  The  event  associated  with  the 
flares  of  August  1972  resulted  In  a  1G3  Increase  In  1  >  4  Ions  at  L  -  2.5 
and  E  ~  0.26  MeV/u.  At  the  same  time,  no  appreciable  Increase  was 
observed  In  the  proton  flux  at  L  ■  2.5.  The  Injected  flux  decayed  with  a 
half-life  of  -  25  days. 


In  addition,  Spjeldnlk  and  Fr1t2  (1981c)  report  that  Z  >  9  nuclei  were 
Injected  by  the  same  flare.  They  found  an  Increase  of  >  1(F  at  L  *  2.5 
from  an  Indetectably  low  flux  prior  to  the  flares.  Assuming  these  nuclei 
were  predominantly  silicon,  they  had  energies  In  the  range  0.4  <  E  <  3.9 
MeV/u.  The  flux  decayed  at  L  «  2.5  with  a  half-life  of  -  30  days.  “This 
decay  time  seemed  to  the  authors,  too  short  to  have  been  controlled  by 
charge  exchange  or  energy  loss  processes. 

The  question  now  Is  what  does  this  Imply  for  the  higher  energies  of 
interest  here.  We  can  only  guess  that  such  injections  must  have  occurred  at 
higher  energies.  If  the  Injections  at  E  i  10  MeV/u  were  as  large  as  those 
reported  at  lower  energies  (and  they  may  well  have  been),  then  heavy  Ions 
were  probably  the  principal  cause  of  soft  errors  arising  from  trapped 
radiations  in  the  entire  Inner  magnetosphere. 

At  higher  energies,  the  decay  times  for  heavy  Ion  Injections  could  be 
many  years,  for  particles  mirroring  near  the  geomagnetic  equatorial  plane. 
With  the  little  we  currently  know  about  these  trapped  energetic  heavy  Ions, 
a  worst  case  model  cannot  be  constructed.  We  simply  have  no  Idea  what  the 
conditions  are  following  a  large  flare  or  how  long  they  persist.  We  can 
only  hope  that  what  has  been  seen  In  the  two  experiments  (Mogro-Campero, 

1972  and  Chan  and  Price,  1975)  Is  typical  and  can  be  extrapolated  to  other 
parts  of  the  magnetosphere  and  other  times. 


With  so  little  experimental  data,  any  model  for  heavy  Ions  In  the 
magnetosphere  Is  highly  speculative.  We  will,  nevertheless,  suggest  a  model 
that  we  hope  will  prove  to  be  sufficiently  conservative. 

For  L  <  3  we  suggest  the  differential  energy  spectrum  for  oxygen  (the 
dashed  line  In  Figure  6.2)  be  adopted,  with  the  ordinate  multiplied  by  72. 
Specifically, 

i  *  £.4  x  10"4  oxygen  nuclei /cluster. sec. MeV/u  E  <  16  MeV/u  (6.5) 


and 

i  -  1.96xl02  E“4,62oxygen  nucle1/cm2ster. sec. MeV/u  E  >  16  MeV/u  (6,6) 

We  take  this  to  be  the  trapped  flux  from  1.2  <  L  <  3.0  for  all  values  of 
B/B0 •  To  obtain  the  spectra  for  other  elements,  Tust“mult1ply  the  oxygen 
spectrum  by  the  measured  relative  abundances  for  the  magnetosphere  In  Table 
2. 


For  L  <  3,  Mogro-Campero  reports  the  ratio  0/C  -  0.5  +  0.4.  This  Is 
more  typical  of  galactic  cosmic  ray  abundances  (see  Table~3). 

No  differential  energy  spectrum  was  measured  In  this  experiment,  hence 
we  can  only  assume  that  the  spectrum  falls  with  increasing  energy  as  the 
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Alfven  criterion  for  stable  trapping  Is  approached.  The  Alfven  criterion 
Is: 


<  1  (6.7) 

where  R  Is  the  magnetic  rigidity  and  11  Is  the  magnetic  field  flux  density. 

This  Is  the  same  factor  that  limits  the  proton  spectrum  at  high 
energies,  so  we  use  the  proton  spectrum  at  L  -  3.0  as  a  model,  from  Sawyer 
and  Vette  (1976): 

«fp  -  E~5.2  (6.8) 

hence  we  propose  to  use 

“  1,6  E'5'2  oxygen  nucle1/cm2ster. sec .MeV/u  for  E  >  13  MeV/u  (6.P) 

for  the  differential  energy  spectrum  of  oxygen.  The  spectra  of  the  other 
elements  are  obtained  from  the  galactic  cosmic  ray  abundances  shown  in 
Table  6.1. 

Besides  nuclei,  there  are  Intense  fluxes  of  electrons  trapped  in  the 
magnetosphere,  especially  In  the  outer  Van  Allen  belt.  In  section  2.2  we 
discussed  the  ways  In  which  electrons  could  cause  soft  upsets  In 
microelectronics  and  concluded  that  we  need  only  be  concerned  about  the  case 
where  the  electron  flux  above  20  MeV  far  exceeds  the  elemental  flux.  There 
Is  a  substantial  flux  of  electrons  above  20  MeV  In  the  Inner  radiation  belt, 
but  the  proton  flux  there  Is  much  greater.  The  electrons  are  the  dominant 
component  In  the  outer  belt,  but  they  are  all  at  energies  below  20  MeV  (see 
Vette  et  al,  1966).  Based  on  these  results,  we  conclude  that  trapped 
electrons  will  not  be  an  Important  cause  of  soft  upsets  on  satellites. 

6.4  Other  Particles  In  the  Magnetosphere 

Besides  stably  trapped  particles,  the  magnetosphere  contains 
quasi -trapped  particles  that  cannot  complete  a  drift  around  the  earth 
without  loss  and  splash  albedo  cosmic  ray  particles.  A  few  experiments  have 
been  carried  out  to  measure  each  of  these  components  (see  Kuznet"'**  ot  al. 
1979,  and  Frledlander  and  Hoppe,  1977).  The  results  give  fluxes  r.nut  are 
very  low  In  comparison  with  other  components  encountered  on  practical 
satellite  orbits,  therefore  they  will  be  Ignored. 
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Fig.  6.2  —  The  orbit-averaged  flux  of  oxygen  nuclei  from  the  Skylab  orbit  aa 
measured  in  a  stack  of  plastic  track  detectors  analyzed  independently  by  Chan 
(3976)  and  Biswas  and  Durgaprasad  (1980).  Also  shown  are  the  galactic  cos¬ 
mic  ray  spectra  (OCR),  geomagnetioally  modulated  for  the  Skylab  orbit  and 
the  anomalous  oxygen  spectrum  similarly  modulated,  assuming  the  oxygen 
nuclei  are  in  the  +1  charge  state. 
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7.0  Cosmic  Ray  Effects  on  Space-Borne  Microelectronics 


In  the  preceding  sections  we  have  developed  a  model  that  describes 
the  particle  environment  through  which  spacecraft  travel.  We  have 
shown  how  one  computes  the  geomagnetic  cutoff  and  uses  It  to  modulate 
the  spectra  of  particles  coming  from  outside  the  magnetosphere. 

Finally,  we  have  suggested  a  model  for  protons  trapped  In  the 
magnetosphere  and  shown  how  this  model  can  be  augmented  to  predict  the 
trapped  fluxes  of  helium  and  heavier  nuclei.  By  following  this 
procedure,  one  can  arrive  at  an  estimate  of  the  differential  energy 
spectra  of  the  various  nuclei  actually  Incident  on  the  spacecraft.  The 
present  section  Is  devoted  to  a  discussion  of  how  these  spectra  are 
used  to  estimate  the  cosmic  ray  effects  on  satellite-borne 
microelectronics  and  what  might  be  done  to  minimize  those  effects. 


7.1  Cosmic  Ray  Transport  Through  the  Spacecraft  Walls 

Since  cosmic  rays  are  Isotropic,  they  may  reach  the  electronics  In 
the  spacecraft  from  all  directions.  Some  of  these  directions  may 
Involve  penetrating  as  little  as  25  mils  of  aluminum,  while  others  may 
Involve  a  few  Inches', 

A  complete  treatment  of  the  cosmic  ray  transport  problem  Involves 
accounting  for  how  the  energy  spectrum  of  each  element  Is  modified  by 
energy  loss  and  nuclear  fragmentation  In  the  spacecraft,  helnrlch 
(1977)  has  carried  out  such  a  calculation  for  the  interplanetary  cosmic 
ray  spectrum.  Ho  shows  that  the  effect  of  shielding  Is  tc  reduce  the 
amplitude  of  the  differential  energy  spectra  and  gradually  shift  the 
peak  of  the  spectrum  to  higher  energies.  It  Is  only  for  minor 
constituent  nuclei,  like  fluorine,  that  fragmentation  can  actually 
Increase  the  flux  above  the  unshielded  level.  Even  then,  this  increase 
Is  only  at  low  energies  and  for  a  limited  range  of  shield  thickness. 

While  the  Interplanetary  cosmic-ray  spectral  shape  is  little 
affected  by  transport  through  shielding,  this  Is  not  true  of  cosmic 
rays  modulated  by  the  geomagnetic  cutoff.  The  effect  of  the  earth's 
magnetic  field  is  to  exclude  the  lowest  energy  cosmic  rays.  After  such 
a  spectrum  passes  through  shielding,  this  low  energy  part  of  the 
spectrum  will  be  restored  by  Ionization  losses  1n£he  shielding. 

Several  workers,  publishing  papers  on  soft  errors,  have  presented 
calculations  of  shielded  cosmic  ray  spectra  which  show  enormous 
increases  In  the  low  energy  cosmic  ray  Intensity  after  propagation 
through  a  modest  shield  thickness.  These  calculations  are  grossly  in 
error!  In  making  such  calculations  care  must  be  taken.  It  is  not 
enough  to  shift  the  cosmic  ray  Intensity  to  a  lower  energy.  The  energy 
Interval  over  which  that  Intensity  is  measured  must  also  be  transformed 
to  the  lower  energy.  When  this  Is  properly  done,  results  similar  to 
those  of  Heinrich  will  be  obtained.  The  only  possibility  for  building 
up  the  spectrum  at  low  energies  is  In  the  case  of  the  hydrogen 
spectrum.  Here  It  Is  possible  for  protons  to  suffer  large  energy 
losses  in  nuclear  collisions.  A  small  number  of  low  energy  protons 
could  be  produced  In  nuclear  reactions  Initiated  by  heavier  ions.  Such 
protons  could  produce  a  low  energy  build  up  In  the  shielded  hydrogen 
cosmic  ray  spectrum. 
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7.2  The  Effects  of  Intense  Ionization 

In  the  previous  paragraphs,  we  saw  how  It  Is  possible  to  obtain  the 
cosmic  ray  spectra  Incident  on  the  microelectronic  devices.  We  will  now 
discuss  the  effects  of  the  Intensely  Ionized  paths  left  by  these  particles. 

Wallmark  and  Marcus  (1962)  were  the  first  to  predict  that  cosmic  rays 
could  have  adverse  effects  on  mlcroel ectronlcs.  These  authors  predicted  two 
effects:  false  signals  resulting  from  cosmic  ray  Ionization  and 
displacement  damage  resulting  In  the  permanent  failure  of  Individual  circuit 
elements.  The  first  of  these  effects  was  observed  In  dynamic  random-access 
memories  (RAMs)  Initially  by  May  and  Woods  (1979).  These  authors  observed 
bits  flipped  by  charged  particles  with  no  permanent  damage  to  the  memory 
cells.  They  called  these  events  "soft  errors".  Following  publication  of 
their  results,  there  has  been  a  flurry  of  papers  describing  soft  errors  In  a 
variety  of  devices. 

May  and  Woods  showed  that  the  soft  errors  they  observed  were  due  to 
alpha  particles  from  uranium  and  thorium  in  the  device  packages.  Clearly, 
whatever  a  shopping  alpha  particle  can  do,  a  stopping  Iron  nucleus  can  do 
much  more  easily)  The  first  published  results  on  tests  with  heavy  Ions  were 
presented  by  Kolaslnskl,  et  a  1 .  (lS‘/9).  These  outhors  reported  trie  results 
of  tests  on  a  variety  of  devices,  many  of  which  are  candidates  for 
spacecraft  electronic  circuits.  Most  of  these  devices  were  found  to  exhibit 
soft  errors.  In  addition  some  of  the  devices  exhibited  latchup  {a  bit  which 
Is  permanently  latched  In  one  state  until  released  by  cycling  the  power 
supplies).  The  mechanism  for  latchup  In  one  of  these  devices,  the  HM-6508, 
has  been  Identified  (Pickel  ar.d  Blandford,  1980b).  It  results  from  "turning 
on"  a  parasitic  circuit  which  Is  operationally  equivalent  to  a 
silicon-controlled  rectifier.  Once  turned  on,  the  circuit  must  be  powered 
down  to  reset  It.  The  "turn  on"  current  comes  from  the  hole-electron  pairs 
generated  by  a  stopping  heavy  Ion. 

There  have  been  quite  a  number  of  efforts  to  model  the  mechanisms  by 
which  cosmic  rays  produce  soft  errors.  May  and  Woods  seem  to  have  correctly 
identified  the  mechanism  In  dynamic  RAMs.  A  soft  error  will  result  whenever 
a  passing  alpha  particle  produces  enough  ionization  to  supply  the  critical 
charge  needed  to  fill  an  empty  potential  well  In  the  transistor-capacitor 
circuit  of  a  dynamic  RAM  memory  call.  Besides  this  mechanism,  there  appear 
to  be  many  others.  For  example,  floating  bit-lines  can  be  discharged,  and 
bistable  flip-flops  may  be  disrupted  by  cosmlc-ray-lnduced  voltage 
trunsients  and  regenerate  In  the  opposite  state.  Pickel  and  Blandford 
(1980a)  have  discussed  models  fo<"  these  and  other  mechanisms  In  detail. 

Besides  soft  errors,  Wallmark  and  Marcus  (I960)  predicted  hard 
failures.  Pickel  and  Blandford  (1980a)  have  reported  hard  failures  due  to 
heavy-ion-inouced  Insulator  punctures  In  MN05  EAROMS.  In  addition,  based  on 
cosmic-ray  track  structure  considerations,  Bradford  11978)  has  predicted 
that  permanent  radiation  damage  will  be  produced  by  cosmic  rays  in  VLSI 
circuits- 
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7.3  Efforts  to  Reduce  Cosmic  Ray  Effects 

There  are  a  number  of  proposals  to  reduce  the  effects  of  cosmic  rays  on 
microelectronics .  The  most  obvious  Is  to  test  available  devices  In  heavy 
ion  beams  and  reject  those  that  are  affected.  This  program  Is  being 
actively  pursued  by  NASA,  SAMSQ,  and  their  contractors  It  Is  becoming 
clear,  from  their  work,  that  certain  circuit  designs  and  device  technologies 
are  much  less  sensitive  than  others. 

A  second  approach  (first  suggested  by  Pickel  and  Blandford,  1978)  .'s  to 
pursue  an  especially  promising  device  technology,  CMOS/SOS.  Because 
hole-electron  pairs  in  the  sapphire  substrates  of  these  devices  have  a  very 
low  mobility,  the  hole-electron  pairs  generated  in  the  substrates  cannot 
readily  diffuse  Into  the  circuits'  potential  wells.  This  makes  it  necessary 
for  cosmic  rays  to  produce  the  required  charge  in  a  much  shorter  path. 

This,  In  turn,  calls  for  higher  stopping  power,  and  thus  heavier  Ions  are 
required  to  produce  soft  errors.  Kolaslnskl  et  al.  (1979)  and  Brucker  et 
al .  (1980)  have  found  some  very  error-resistant  devices  of  this  type. 

Another  approach,  tried  with  some  success,  is  to  adjust  the  device 
parameters.  Peeples  and  Every  (1980)  show  that  soft  error  rates  can  be 
reduced  by  increasing  supply  voltages.  Increasing  voltages,  however, 
increases  the  device  hard-failure  rate  and  is  therefore  of  limited  ur,a  In 
controlling  soft  errors. 

Throughout  the  test  results  published  to  date  there  Is  a  disturbing 
correlation.  The  more  dense  and  lower  power  a  device  Is  made,  the  more 
sensitive  It  becomes  to  cosmic  rays.  Tne  most  advanced  devices,  64K  dynamic 
RAMs,  256K  CCD's  and  VLSI  circuits  are  highly  desirable  candidates  for 
satellite  applications,  but  they  are  also  among  the  most  sensitive  devices 
to  cosmic  rays.  A  potential  solution  to  the  problem  Is  a  fault-tolerant 
space-borne  computer.  Designs  for  such  computers  have  recently  reported  by 
Retzler  (1980)  and  Masson  (1980).  While  these  designs  appear  capable  of 
tolerating  moderate  error  rates,  they  Impose  a  large  penalty  In  weight, 
power  and  cost. 

It  seems  clear  that  regardless  of  technique  employed  to  reduce  soft 
errors,  an  accurate  means  must  be  developed  to  estimate  soft  error  rates 
since  large  design-penalty  trade-offs  will  hinge  on  this  parameter. 

7.4  Shielding  Against  Cosmic  Rays 

Spacecraft  shielding  is  not  generally  effective  against  cosmic  rays. 
There  are,  however,  specific  instances  in  which  It  can  clearly  reduce  or 
1 ncrease  the  soft  error  rate. 

Suppose  a  spacecraft  Is  In  a  13,000  km  circular  orbit  with  a  low 
Inclination.  This  spacecraft  would  sample  the  particles  trapped  between 

1,  *  3  and  L  =  3.5  while  always  having  a  geomagnetic  cutoff  (from  eq.  5.5)  of 
p  ~  I6/L2  -  1.3  to  1.8  GV 

i.e.  180  to  320  MeV/u.  Wow  suppose  the  electronics  are  lightly  shielded  on 
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one  side,  say  25  mils  of  aluminum,  and  the  components  are  sensitive  to 
helium,  or  even  carbon  nuclei.  Under  these  conditions,  the  trapped  helium 
and  CNO  (carbon,  nitrogen  and  oxygen)  fluxes  at  low  energies,  tor  these 
altitudes,  could  be  a  major  cause  of  soft  errors.  There  are  two  reasons: 
first,  helium  fluxes  at  this  altitude  are  comparable  to  S1(p,a)-produced 
secondary  helium  fluxes  from  trapped  protons  (see  Peterson,  1 980 ) ;  second, 
provided  the  shielding  of  the  electronics  does  not  exceed  about  0.6  Inches 
In  any  direction,  cosmic  ray  CNO  pass  through  the  electronics  too  fast  to 
have  the  required  stopping  power.  In  this  case,  only  trapped  CNO  can 
contribute.  Now  by  Increasing  the  shielding  to,  at  least,  100  mils  In  all 
directions,  the  trapped  He,  C,  N,  and  0  nuclei  will  probably  be  excluded  and 
the  error  rate  reduced. 

This  example  Is  a  rather  restricted  case.  A  much  more  general  case  Is 
for  a  spacecraft  In  Interstellar  space,  or  nearly  so.  Suppose  that  the 
electronics  are  sensitive  only  to  stopping  iron,  and  are  required  to  have 
low  soft  error  rates,  even  during  an  anomalously  large  solar  flare.  If  the 
electronics  can  be  made  to  tolerate  such  a  flare,  galactic  cosmic  rays  will 
never  produce  an  unacceptable  soft  error  rate.  The  spectra  of  all  flares 
fall  steeply  from  the  lowest  energies,  so  that  the  addition  of  shielding  is 
especially  effective  In  reducing  the  flux  at  the  electronics.  Furthermore, 
since  Iron  has  a  relatively  short  range  [(£6/262)  X  the  range  of  a  proton 
at  the  same  velocity]  the  shielding  will  be  crucially  effective  In  reducing 
the  iron  flux  at  she  electronics.  Clearly,  1;.  .uch  a  case,  shielding  can  be 
beneflc  i  ,il . 

Sh'  Iding  does  not  always  help.  Consider  a  low  altitude,  400  km 
circular  orbit  with  a  low  Inclination  such  that  the  spacecraft  never  goes 
beyond  L  *  3.  Further,  suppose  that  the  electronics  are  only  sensitive  to 
stopping  Iron.  Because  o*  this,  the  trapped  protons  encountered  In  the 
South  Atlantic  anomaly  cannot  produce  soft  errors,  because  the  worst  they 
can  do  Is  to  produce  recoil  silicon.  There  are,  of  course,  heavier  nuclei 
trapped  in  the  radiation  belts,  even  iron,  but  these  are  probably  rare  and 
low  in  energy.  A  minimum  of  100  mils  of  aluminum  will  probably  stop  all  of 
these  rare  nuclei.  Such  a  satellite  would  not  be  troubled  by  soft  errors 
provided  it  Is  not  too  well  shielded.  The  lowest  geomagnetic  cutoff 
encountered  on  the  orbit  is  320  MeV/u  (from  eq.  5.5).  At  this  energy  an 
iron  nucleus  has  a  range  of  880  mils  in  aluminum.  If  there  is  a  direction, 
any  direction,  looking  outward  from  the  electronics  through  that  much 
material,  a  cosmic  ray  Iron  arriving  along  such  a  direction  could  stop  In 
the  electronics  and  cause  an  error.  All  It  takes  is  a  structural  beam 
"aimed"  at  the  electronics.  Because  nearby  objects  subtend  greater  solid 
angles,  the  electronics  may  shield  themselves.  The  habit  many  designers 
have  of  placing  components  in  a  neat  little  row  can  build  up  a  large 
shielding  thickness  along  that  row,  especially  if  It.'s  a  row  or,  say, 
tantalum  capacitors! 

From  these  examples  the  reader  can  see  that,  while  shielding  can  be 
important,  Its  effects  depend  on  specific  mission  requirements,  orbits  and 
device  sensitivities.  Indeed,  estimating  soft  error  rates  Is  not  a  simple 
business! 
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l.o  The  Relative  Importance  of  the  Various  Components  in  the  Near-Earth 
Particle  Environment. 

From  what  has  been  said  up  to  now,  It  Is  probably  clear  that  given  the 
right  set  of  circumstances  (l.e.  orbit,  tolerable  error  rate,  shielding,  and 
device  sensitivity),  any  of  the  components  we  have  described  could  be  the 
dominant  source  of  soft  errors.  In  the  real  world,  however,  there  are 
"popular  orbits,"  actual  devices  and  typical  missions.  In  the  following 
paragraphs  we  will  discuss  how  the  environment  affects  some  of  these 

Consider  a  28.5°  Inclination  circular  orbit  at  300-500  km.  This  Is  a 
typical  Eastern  Test  Range  minimum-energy  orbit  for  shuttle  sortie  or 
shuttle-deployed  free  filer  payloads.  For  electronics  sensitive  to  eny 
secondary  particle  that  a  proton  can  produce  (see  Peterson,  1980).  the  South 
Atlantic  anomaly  will  be  the  problem  area,  and  shielding  will  not  help  much 
against  the  protons.  In  this  part  of  the  magnetosphere,  th£  proton  spectrum 
Is  very  hard  (see  Sawyer  and  Vette,  1976).  A  little  shielding  would  help 
«ga1r,st  heavier  trapped  nuclei,  but  their  flux  In  this  orbit  may  be  too  low 
to  matter  anyway. 

If  the  electronics  In  this  payload  cannot  be  upset  by  proton-produced 
secondaries  then  something  like  stopping  iron  nuclei  will  be  necessary  and 
the  soft  error  rate  will  be  dramatically  lower.  The  cosmic  ray  cutoff  for 
this  orbit  Is  ~  1000  KeV/u.  At  that  energy,  It  takes  5  Inches  of  aluminum 
to  stop  a  cosmic  ray  Iron  nucleus.  Even  then  only  about  20  per  cant  of  them 
will  survive  to  come  to  rest  without  fragmenting  in  a  collision  with  *n 
aluminum  nucleus.  Assuming  the  electronics  are  shielded  with  much  less  than 
5  Inches  of  aluminum  or  equivalent,  the  only  things  left  to  worry  about  ore 
trapped  heavy  nuclei,  especially  Iron.  They  can  probably  be  dealt  with  by 
100  mils  or  so  of  shielding.  All  that  remains  then  is  cosmic  ray  re-entrant 
albedo,  l.e.  cosmic  rays  that  skip  off  the  atmosphere,  thus  evading  the 
geomagnetic  cutoff  and  arriving  at  the  spacecraft  with  a  lower  energy  than 
1000  NieV/u.  These  events  are  very  rare  and  can  probably  be  neglecteci. 

There  Is  always  the  anomalous  component,  If  it's  singly  ionized  and  if 
besides  oxygen.  It  contains  Ions  like  Iron.  It  couid  be  a  problem,  but  the 
charge  state  of  the  anomalous  component  is  unknown.  For  the  present  we 
don't  know  If  It  could  contribute,  so  we  can  be  optimistic  and  ignore  It. 

Solar  flares  could,  In  principle,  cause  soft  errors  in  this  payload,  but 
they  are  subject  to  almost  the  same  geomagnetic  cutoff  as  cosmic  rays.  I* 
the  electronics  are  so  lightly  shleleed  that  cosmic  rays  don't  affect  them, 
flares  wll 1  not  either. 

Next,  consider  a  low  altitude  circular  orbit  of  300  to  500  km  again,  but 
this  time  Inclined  at  90o,  l.e.  polar.  For  electronic  components  that  are 
sensitive  to  trapped  protons,  the  problems  In  the  South  Atlantic  anomaly 
will  be  worse.  Instead  of  only  brushing  through  the  northern  side  of  the 
anomaly,  as  In  the  previous  example,  this  payload  will  pass  through  the 
heart  of  the  anomaly.  The  error  rate  experienced  In  the  anomaly  will  depend 
critically  on  altitude.  The  proton  flux  considered  here,  scales 
approximately  as  altitude  to  the  5th  power!  At  a  300  km  altitude,  It  Is 
almost  possible  to  sneak  through  under  the  anomaly. 

The  southern  port  of  the  anomaly,  at  L  values  between  2  and  2.5, 
probably  contain  enough  energetic  helium  to  contribute  directly  to  the  error 


68 


rate.  This  1 s  south  of  the  most  Intense  part  and  the  error  rate  will  be 
lower  here,  whatever  the  cause.  Extra  shielding.  100  mils  or  so  could  help 
here,  but  It's  probably  not  worth  the  weight  and  bother. 

While  It  may  be  possible  to  sneak  through  under  the  anomaly,  there  is  no 
way  to  "sneak  through"  under  the  polar  clefts.  When  the  spacecraft  passes 
over  the  magnetic  poles,  the  geomagnetic  cutoff  goes  practically  to  zero. 
Depending  on  the  orbital  altitude  and  the  sensitivity  of  the  electronics, 
cosmic  rays  coming  In  over  the  poles  may  cause  a  higher  error  rate  than  the 
South  Atlantic  anomaly,  for  electronics  not  sensitive  to  proton-produced 
secondaries,  cosmic  rays  will  surely  be  the  major  cause  of  soft  errors, 
regardless  of  altitude. 

If  this  payload  Is  designed  to  operate  during  a  major  flare,  Its 
electronics  will  easily  tolerate  cosmic  rays  and  the  South  Atlantic 
anomaly.  The  payload  should  be  as  heavily  shielded  as  possible,  If  the 
electronics  are  even  sensitive  to  very  heavy  nuclei.  As  noted  earlier,  the 
shielding  Is  especially  effective  In  stopping  these  nuclei  and  reducing 
their  Intensity  In  the  steeply  falling  spectrum  of  a  solar  flare. 

As  a  final  example,  we  consider  a  payload  in  synchronous  orbit  with  0° 
Inclination.  The  geomagnetic  cutoff  at  synchronous  orbit  Is 

P  -  16/(6. 6)2  0.364  C V 

or  -  15.2  MeV/u.  At  this  energy,  protons  and  helium  nuclei  have  a  range 
of  50  mils  of  aluminum;  for  heavier  nuclei  It  will  be  even  shorter. 

Clearly,  there  is  no  realistic  way  to  shield  the  electronics  so  lightly  that 
cosmic  ray  heavy  ions  do  not  stop  In  the  circuits.  For  cosmic  rays  and 
solar  flares,  It's  just  as  though  the  spacecraft  were  In  1 nterplanetary 
space. 

As  Heinrich  (1977)  has  shown,  ever  10  g/cm2  (■  1.5  inches  of  aluminum) 
doesn't  do  much  to  alter  the  cosmic  ray  spectra  above  10  MeV/u.  Clearly, 
shielding  effectively  aga  nst  galactic  cosmic  rays  Is  hopeless.  Some 
researchers  (Oernart  and  Stekly,  1964)  have  suggested  using  a  strong 
magnetic  field  coupled  with  shielding  to  protect  against  cosmic  rays.  We 
judge  such  a  proposal  as  heavy,  expensive  and,  worst  of  all,  dubious,  since 
such  an  artlflcal  magnetic  field  In  space  would  develop  Its  own 
magnetosphere  filial!  with  trapped  particles.  The  cure  could  be  worse  than 
the  disease! 

Since  this  satellite  Is  unprotected  by  the  magnetosphere,  It  Is  exposed 
to  solar  flare  particles.  A  decision  has  to  be  made  here.  It's  much  easier 
to  design  for  a  tolerable  error  rate  from  cosmic  rayc  than  from  a  major 
solar  flare.  If  the  mission  can  tolerate  being  shut  down  during  a  major 
flare,  the  design  problem  will  be  much  easier.  This  may  represent  a  loss  of 
only  ~  2  per  cent  of  the  mission  time  and  this  has  to  be  traded  against 
the  80-90  per  cent  da'sa  recovery  rates  that  are  typical  of  many  space 
missions.  In  short,  losing  data  during  major  flares  doesn't  cost  much  in 
mission  time. 
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There  are  some  payloads  that  must  continue  to  operate  correctly  during 
flares,  unlnterruptable  military  missions  and  scientific  experiments  to 
study  flares,  for  Instance.  For  these  payloads  one  must  reckon  with  the 
model  In  section  4.0.  The  mean,  worst  case,  and  anomalous  case  flare 
spectra  are  shown  In  figures  4.2  and  4.3.  All  these  spectra,  unlike 
galactic  cosmic  rays,  are  monotonlcally  and  rapidly  decreasing  functions  of 
energy.  By  comparing  the  spectra  In  figure  4.3  with  the  cosmic  ray  proton 
spectrum,  figure  2.2,  we  see  that  the  peak  flare  flux  falls  below  the  cosmic 
ray  background  between  200  MeV/u  and  -  1000  MeV/u  depending  on  flare 
size.  Since  flares  are  not  generally  as  rich  in  the  heavier  elements  as 
cosmic  rays,  this  crossover  may  occur  at  somewhat  lower  energies  when 
comparing  iron  spectra. 

Now  let  us  assume  that  the  spacecraft  designer,  faced  with  the  problem 
of  controlling  the  soft  error  rate  during  a  flare,  resorts  to  circuits  that 
are  relatively  "cosmic-ray  hard"  l.e.  only  very  heavy  stopping  nuclei  are 
capable  of  causing  errors.  He  may  be  able  to  use  components  that  are  so 
insensitive  to  soft  errors  that  the  problem  can  he  controlled  In  that  way 
alone.  This  solution,  however,  may  not  be  possible  within  the  weight,  power 
and  cost  constraints.  If  components  which  are  just  sensitive  to  stopping 
very  heavy  nuclei  must  be  used,  then  shielding  can  help.  Suppose  an 
anomalously  large  flare  occurs  with  a  peak  spectrum  of  the  kind  shown  in 
figure  4.3.  The  stopping  power  of  an  Iron  nucleus  reaches  Its  peak  at  about 
2  MeV/u.  Under  25  mils  of  aluminum,  the  flux  of  Iron  at  this  energy  will  be 
t.l  x  105  particles/m2  ster.  sec.  MeV/u.  If  100  mils  of  shielding  Is 
provided,  the  flux  behind  the  shield,  again  at  2  MeV/u  will  be  E.8  x  104 
part1cles/m2  ster.  sec.  MeV/u.  This  Is  an  order  of  magnitude 
improvement.  Additional  shielding  will  further  reduce  the  flux  at  2  MeV/u, 
but  not  as  rapidly.  It  would  be  necessary  to  add  several  Inches  of  aluminum 
to  reduce  the  flare  flux  to  near  cosmic  ray  background  levels.  While  this 
goal  Is  unreasonable,  some  trade-off  on  shielding  should  be  possible. 

In  these  examples  we  have  tried  to  show  that  galactic  cosmic  rays,  solar 
flare  particles,  trapped  protons  and  possibly  even  trapped  heavy  nuclei  can 
be  the  dominant  causes  of  soft  errors  under  the  right  conditions.  Unless 
the  mission  objectives  can  be  economically  accomplished  with  components  tnat 
are  practically  "Immune"  to  soft  errors,  it  appears  that  an  accurate 
estimate  of  soft  error  rates  will  be  a  critical  parameter  In  deciding  the 
various  tradeoffs  that  must  be  made  In  mission  planning. 
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6.0  Conclusions  and  Recommendations 


The  published  data  and  their  interpretation  described  in  section  7.2 
make  it  clear  that  Intensely  ionizing  particles  can  produce  soft  errors  and 
related  phenomena.  Furthermore  it  is  well  established  (May  and  Noods,  1979) 
that  there  is  a  threshold;  a  critical  charge  that  must  be  liberated  by  the 
ionizing  particle  in  a  very  short  path  length.  Below  this  threshold,  soft 
errors  do  not  occur.  For  fixed  device  dimensions  and  characteristics  this 
threshold  becomes  a  threshold  on  stopping  power,  which  varies  approximately 
as  (the  particle's  charge)2/(the  particle's  velocity)2.  Because 
electrons  are  picked  up  from  the  medium  by  a  stopping  ion,  there  is  a 
maximum  stopping  power  reached  by  each  ion  just  before  coming  to  rest.  If 
the  threshold  is  higher  than  the  maximum  stopping  power  of  an  ion,  then  it 
cannot  produce  a  soft  error;  a  more  highly  charged,  i.e.  heavier  ion  will  be 
required. 

The  critical  charge  is  a  key  device  parameter  that  helps  determine  a 
device's  sensitivity  to  soft  errors.  For  satellite  applications,  an  effort 
must  be  made  to  select  devices  that  are  insensitive,  I.e.,  have  a  high 
critical  charge.  This  means  in  turn,  that  only  heavy  or  very  heavy  ions  can 
produce  errors  in  such  devices. 

6.1  Galactic  Cosmic  Rays 

We  have  made  a  comprehensive  search  of  the  published  data  on  cosmic  rays 
and  developed  an  analytic  model  for  the  cosmic  ray  environment.  This  model 
describes  the  differential  energy  spectra  for  cosmic  rays  at  all  energies 
greater  than  10  MeV/u  and  for  all  nuclei  in  the  periodic  table  up  to 
nickel.  Me  judge  that  the  model  fits  the  data  to  +  20  per  cent  for 
hydrogen,  helium  and  the  more  abundant  elements  up  to  neon  at  energies  below 
104  MeV/u.  For  higher  energies  and  less  abundant  elements  the  fits  are 
probably  accurate  to  +  50  per  cent.  In  one  Important  case,  the  available 
data  were  not  adequate  to  define  the  model  well.  This  is  the  differential 
energy  spectrum  of  iron,  figure  2.7.  There  are  no  data  on  the  spectrum  at 
solar  maximum  for  energies  below  900  Mev/u.  Secondly  the  data  below  200 
Mev/u  are  very  sparse,  and  only  two  experiments  give  results  below  100 
Mev/u.  To  emphasize  the  Importance  of  the  iron  spectrum  below  100  Mev/u,  we 
note  that  iron  nuclei  with  energies  between  50  and  100  Mev/u  are  just 
stopping  as  they  exit  aluminum  shielding  with  thicknesses  in  the  range  of  25 
to  120  mils.  It  is  Just  these  stopping  iron  nuclei  that  are  the  most 
effective  in  producing  soft  errors. 

Except  for  the  data  on  the  iron  spectrum,  the  galactic  cosmic  ray  data 
base  we  used  is  reasonably  complete.  The  errors  In  the  analytic  model  we 
suggest,  are  generally  smaller  than  some  of  the  other  sources  of  error  in 
the  problem. 

The  model  Includes  the  effects  of  solar  modulation,  by  sinusoidally 
interpolating  between  the  solar  maximum  and  solar  minimum  analytic-model 
spectra.  When  flux  levels  must  be  predicted  in  the  future,  such  a  model  has 
to  be  employed.  If  we  take  our  model  for  solar  cycle  variations,  and 
compare  it  with  the  past,  Figure  2.14,  we  see  that  it  is,  at  best,  only  a 
crude  fit  to  the  data.  From  this  figure,  it  Is  clear  that  a  major  source  of 
error  in  these  predictions  is  our  inability  to  predict  future  levels  of 
solar  modulation. 
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The  prediction  of  future  levels  of  solar  modulation  has  received  very 
little  attention  by  the  scientific  community  and  yet  It  appears  to  be  the 
major  source  of  error  In  our  model  of  galactic  cosmic  rays.  Comparing  the 
analytic  fit  to  the  solar  maximum  and  minimum  helium  spectra  at  50  Mev/u, 
figure  2.3  with  the  solar  modulation  fit  and  data,  figure  ?!.14,  we  see  that 
our  predictions  could  be  wrong  by  a  factor  of  5  at  solar  maximum,  and  a 
factor  of  2.6  at  solar  minimum.  The  overall  accuracy  of  the  model 
predictions  Is  probably  no  better  than  +  100  per  cent  at  low  energies.  Also 
It  should  be  noted  that  these  low  energy  cosmic  rays  are  just  stopping  as 
they  enter  the  electronic  components,  l.e.  the  ones  most  effective  In 
causing  soft  errors. 

8.2  The  Anomalous  Component 

This  1  Ittle-understood  feature  1  ri  the  cosmic  ray  spectra  of  certain 
elements  was  discussed  In  Section  3.0.  its  presence  In  the  Interplanetary 
medium  near  earth  was  obvious  from  1972  until  1978.  Whether  this  feature 
will  reappear  during  the  1983-1989  time  frame  Is  unknown.  Indeed,  since  It 
never  went  away  In  the  outer  solar  system,  It  may  always  be  with  us  at  some 
level  here,  near  earth.  The  greatest  mystery  about  the  anomalous  component 
Is  whether  It  Is  singly  Ionized.  If  It  Is  singly  ionized,  Its  presence 
Inside  the  magnetosphere  would  be  much  more  widespread.  Since  the  anomalous 
component  Is  a  feature  of  the  low  energy  part  of  the  cosmic  ray  spectrum, 
these  nuclei  will  be  most  effective  in  causing  soft  errors.  Tnelr  wide 
spread  access  to  the  Inner  magnetosphere  could  change  the  soft  error  rates 
on  some  satellites  dramatically .  Clearly  the  charge  state  of  the  anomalous 
component  should  be  established  as  soon  as  possible  because  of  Its  Impact  on 
the  predictions  of  soft  error  rates. 

8.3  Solar  Flares 


About  2  per  cent  of  the  time,  a  major  solar  flare  Is  In  progress,  with 
an  associated  solar  particle  event  In  the  Interplanetary  medium.  When  these 
events  occur,  they  will  be  the  dominant  source  of  soft  errors  for  satellites 
that  are  not  well  shielded  by  the  earth's  magnetic  field.  Not  all 
satellites  will  have  to  operate  during  such  events,  but  many  will  have  to 
continue  reliable  operation  under  the  heaviest  of  solar  weather  conditions. 

Section  4.0  addressed  the  problems  to  be  faced  In  planning  for  solar 
flares.  The  fundamental  problem  with  flares  Is  the  uncertainty  wc  have 
about  their  frequency,  size  distribution,  and  composition  variability. 

There  Is  a  substantial  research  effort  In  the  field  of  solar  physics  being 
supported  by  NASA.  Several  of  these  research  programs  are  aimed  at 
energetic  solar  particle  measurements,  but  mostly  at  energies  too  low  to 
affect  satellite  electronics  directly  .  The  Office  of  Naval  Research  Is 
sponsoring  two  solar  particle  experiments  from  the  University  of  Chicago, 
but  we  are  unaware  of  any  other  research  In  this  area  supported  by  the  DoD. 

Solar  flares  have  been  recognized  as  the  cause  of  other  problems  for 
satellites,  In  particular,  and  communications  In  general,  and  are  the 
subject  of  ongoing  research.  The  unanswered  questions  of  special  relevance 
to  the  soft  error  problem  are:  1)  how  can  we  best  predict  the  size  and 
frequency  of  flares,  2)  how  can  we  best  characterize  their  energy  spectra, 
and,  3)  perhaps  most  Important,  how  can  we  best  describe  the  solar  flare 
composition  and  Its  variability. 
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Besides  tne  major  solar  flares,  there  are  small  flares  and  particle  flux 
enhancements  occurring  much  more  frequently.  The  uncertainty  this  Introduces 
Into  the  flux  levels  at  low  energies  was  discussed  In  section  2.1  based  on 
unpublished  IMP  satellite  data  supplied  us  by  Dr.  Robert  Pyle  of  the 
University  of  Chicago.  There  have  been  no  publications,  to  our  knowledge, 
discussing  the  variability  of  low  energy  fluxes.  The  data  base  for  such  a 
study  has,  no  doubt,  already  been  collected  by  a  number  of  groups  doing 
research  In  these  areas.  All  that  remains  Is  to  analyze  the  existing  data. 

8.4  Trapped  Particles 

There  exists  a  standard  model  for  the  trapped  proton  environment  (Sawyer 
and  Vette,  1970).  At  present,  this  model  describes  a  static  magnetosphere 
with  two  states  corresponding  to  solar  maximum  and  minimum  conditions. 
Obviously,  the  magnetosphere  Is  dynamic  with  both  short  term  and  long  term 
variations  in  the  trapped  particle  population.  No  model  has  been 
constructed  to  describe  these  variations,  primarily  because  the  data  base 
does  not  exist.  In  spite  of  this  deficiency,  studies  of  the  trapped 
population  have  been  very  limited  In  recent  years.  Unless  this  trend  Is 
reversed,  It  Is  doubtful  that  a  date  base,  suitable  for  such  a  model,  will 
ever  be  collected. 

Besides  protons,  there  are  other  nuclei  trapped  In  the  magnetosphere. 
Helium,  carbon,  and  oxygen  nuclei  have  all  been  Identified  and,  no  doubt  all 
heavier  elements  In  the  periodic  table  appear  to  be  present  as  well.  While 
there  have  been  a  number  of  experiments  studying  these  nuclei  at  very  low 
energies,  only  four  experiments  have  been  performed  to  examine  the  fluxes  at 
energies  above  10  MeV/u.  It  Is  our  opinion  that  the  results  of  these 
experiments  are  In  disagreement  with  a  steady-state  radial  diffusion  model 
that  has  the  solar  wind  as  a  source,  though,  this  model  seems  to  explain 
most  results  reasonably  well  at  lower  energies.  From  the  results  of  these 
four  experiments,  It  appears  that  trapped  heavy  nuclei  may  contribute  to  the 
soft  error  problem,  at  least,  In  some  parts  of  the  magnetosphere  though 
possibly  only  for  lightly  shielded  satellites. 

In  light  of  the  mysterious  origin  of  these  nuclei,  the  exceedingly 
scanty  experimental  data,  and  their  potential  Impact  on  spacecraft  design, 

It  will  be  necessary  to  carry  out  additional  experiments  to  confirm  the 
earlier  results  and  extend  the  measurements  to  obtain  a  firmer  basis  for 
accessing  the  contribution  of  trapped  heavy  Ions  to  the  soft  errors  observed 
on  satellites. 

Some  final  words  of  caution  need  to  be  mentioned  about  the  particle 
environment.  In  general.  We  have  paid  very  little  attention  to  nuclei 
heavier  than  Iron,  yet  cosmic  rays  have  been  detected  with  charges  up  to 
-  96  (the  element  Curium  or  thereabouts).  While  these  nuclei  are  rare, 
they  do  exl st  In  cosmic  rays  and  they  are  probably  capable  of  upsetting 
mlcroeTect'ronlc  circuits  normally  thought  of  as  immune  to  soft  errors. 

Claims  of  Immunity  to  cosmic  ray  effects  for  electronic  circuits  must, 
therefore  be  regarded  with  some  skepticism. 

As  we  explained  In  the  introduction,  when  the  data  base  Is  Inadequate  to 
describe  the  environment,  a  credible  worst-case  model  might  be  so 
conservative  that  It  would  result  In  undue  preparation  for  conditions  that 
never  occur.  We  have  chosen,  instead,  to  be  optimistic  and  present  the 
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mildest  environment  that  the  data  will  support.  The  user  should  then  be 
reasonably  guaranteed  of  actually  experiencing  an  environment  close  to  that 
predicted.  He  can  be  confident  that  his  preparations  and  precautions  were 
not  wasted  effort.  It  Is,  of  course,  possible  that  actual  conditions  could 
be  far  more  severe  than  predicted.  Until  experiments  actually  show  that 
speculations  concerning  these  conditions  are  correct,  spacecraft  designers 
must  simply  take  a  chance. 
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APPENDIX  1 

THE  ANALYTIC  MODEL  FOR  THE  CHARGED  PARTICLE  ENVIRONMENT 

In  the  paragraphs  that  follow  we  will  present  a  recipe  for  the 
near-earth  charged  particle  environment.  The  recipe  consists  of  a  set  of 
equations  and  tables  for  computing  the  differential  energy  spectra  of  the 
most  Important  charged  particle  populations  in  the  earth's  vicinity.  These 
equations  were  devised  to  fit  the  data  and  are  intended  to  have  no  physical 
Interpretation.  This  analytic  recipe  may  be  easily  programmed  for  a  digital 
computer  of  almost  ar\y  size,  and  Is  Intended  to  become  a  subroutine  In  a 
program  which  will  be  used  to  estimate  the  soft  error  rates  In  satellite- 
borne  electronics. 

Galactic  cosmic  rays  consist  of  electrons  and  the  nuclei  of  all  the 
elements  in  the  periodic  table,  the  first  28  elements  are  the  most  Important 
for  cosmic  ray  effects  on  microelectronics.  These  particles  are  f»-om 
outside  the  solar  system  and  their  flux  at  low  energies  Is  anti -correlated 
with  solar  activity  (i.e.  more  cosmic  rays  at  solar  minimum).  The 
differential  energy  spectra  In  particles  per  square  meter  -  steradlan  r 
second  -  million  electron  volts  per  atomic  mass  unit  (I.e.  particles/nr 
ster.sec.  MeV/u)  are  given  In  the  following  paragraphs. 

The  spectra  for  protons  (hydrogen  nuclei),  a-particles  (helium  nuclei), 
and  iron  nuclei  are  given  below: 

F(E,t)  *  A(E)sin[W(t-t0)]  +  B(E)  (1) 

where 

W  =  0.576  radian/years, 

t0  =  1950.6  A.D.  date, 
and,  t  =  current  date  in  years, 

B(E)  =  0.5  [fmi*n  (E)  +  f^ax  (E)],  (2) 

and 

A(E)  =  0.5  [fmin  (E)  "  ^max  ^) 

fmax  and  fmi-n  differ  only  by  the  choice  of  constants  in  the  equation, 

f(E)  =  10m  (E/E0)a  (4) 

where  E  is  in  MeV/u 
where, 

a=  a0  |  1  -  expC-Xj  (log10E)b] j  (5) 

and 

m  =  Cj_  exp[-X2  (log10E}2]  -C2  (6) 
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The  values  of  the  constants  a0l  Eq,  b,  Xi,  X?»  Cl  and  Iron 

dven  in  Table  1  for  each  of  the  elements  hydrogen  (H),  helium  (He)  and  iron 

(Fo)  for  the  conditions  of  solar  maximum  and  solar  minimum, 

TABLE  1.  Constants  used  Eqs.  (1-6)  to  Compute 
Energy  Spectra  of  H,  He  and  Fe  at  Solar  Maximum  and  at  Solar  Minimum 


Element 


H-mln 
H-max 
He -min 
He -max 
Fe-mln 
Fe-max 


-2.2 

-2.26 

-2.25 

-2.70 

-2.70 


1.175x10s 

1.175x10s 

7.94xl04 

7.94X104 

1.1x10s 

1.1x10s 


The  Differential  energy  spectra  for ,ca^"  (SI)1™ 

M  %  sir.*1*4™ 

[obtained  from  eq.  (1)]  by  the  appropriate  scaling  factor  In  Table  2. 

TABLE  2:  The  Ratio  of  the  Abundance  of  Various  Nuclei  to  Helium 


Element 


Ratio 

2.5xl0"2 

2.3xl0"2 

4.1xl0’4 

3.5xl0"3 

7.0xl0"4 


Element 


Ratio 

4.7xl0“3 
8. 3xl0"4 
3. 5xl0"3 
2.0xl0“4 
7.4xl0“4 


The  differential  energy  spectra  for  calcium  (Ca),  cobalt  (Co)  and 
nickel  (Ni )  are  obtained  by  multiplying  the  Iron  spectrum  [obtained  from  eq 
(1)]  by  the  scale  factors  listed  in  Table  3. 


TABLE  3:  The  Ratios  of  the  Abundance  of 
Various  Elements  to  Fe 


Element 

Ratio 

Ca 

2.3X10"1 

Co 

6x10” 3 

Ni 

4.8xl0”2 

The  spectra  of  the  elements  lithium  (Li),  beryllium  (Be)  and  boron  (B) 
are  obtained  from  the  helium  spectrum,  F^e.  modified  by  the  equation: 


F*= 


0.0142  FHe, 


0.67  E-0.443  FHe, 


E  <  6000  MeV/u 


(7) 

E  >  6000  MeV/u 


to  obtain  the  combined  spectrum  of  (Li  +  Be  +  B).  Eq.  (7)  is  then  multiplied 
by  the  ratios  in  Table  4  to  obtain  the  individual  elemental  spectra. 


TABLE  4:  The  Relative  Fractions  of  Li,  Be,  and  B  in  the 
Combined  Total  Abundance  Li  +Be  +  B 


Element  Ratio 

Li  0.33 
Be  0.17 
B  0.5 


The  spectrum  of  the  element  nitrogen  (N)  is  obtained  by  modifying  the 
helium  spectrum,  F^e  as  shown  below: 

Ffj  =  {  6.4  xlO'3exp[-.4(LogioE  -  3.15)2] 

+  5.6  x  10-3exp[-.9(LogioE  -  0.8)2]}  F^e  (8) 

where  E  is  in  MeV/u. 


The  spectra  for  the  elements  chlorine,  (Cl),  argon  (Ar),  potassium  (K), 
scandium  (Sc),  titanium  (Ti),  vanadium  (V),  chromium  (Cr)  and  manganese  (Mn) 
are  all  obtained  by  modifying  the  iron  spectrum  Fpe  as  shown  below: 


F*=  Q(E)  FFe 

Q(E)  =  16[l-exp( -. 126  E0*4)]r-0.33 


where  E  is  in  MeV/u. 
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(9) 

(10) 


Finally  F*  the  sub-iron  spectrum  [eq.  (9)]  Is  multiplied  by  the  appropriate 
ratio  In  Table  5  to  obtain  the  Individual  elemental  spectra. 

TABLE  5.  The  Fractional  Abundance  of  each  Element 
In  the  sub-iron  group 


Element 

Ratio 

Element 

Ratio 

Cl 

.07 

Ti 

.14 

Ar 

.13 

V 

.07 

K 

.09 

Cr 

.14 

Sc 

.05 

Mn 

.10 

The  recipe  given  above  is  correct  for  quiet  periods  In  the 
interplanetary  medium  when  only  the  galactic  cosmic  rays  are  present.  These 
conditions  are  often  disturbed,  especially  at  low  energies,  by  small  solar 
flares,  co-rotating  events,  etc.  To  allow  for  typical  disturbed  conditions, 
we  recommend  that,  below  100  MeY/u,  a  worst-case  spectrum  be  employed.  With 
90  per  cent  confidence,  the  particle  flux  should  never  be  more  intense  than 
described  by  this  case. 

To  construct  the  worst-case  spectrum  for  protons,  compute  the  "H-min" 
spectrum  (using  Eq.  4)  and  then  compute  Fn-worst  as  shown  below: 

^H-worst  «  [1897e"E/9-66  +  i.64]FH-nrfn  (H> 

This  applies  for  E  £  100  MeV.  For  higher  energies,  you  may  use  the  galactic 
cosmic  ray  spectrum  for  the  appropriate  mission  time  t  [in  Eq.  (1)]. 

In  like  manner,  the  solar  minimum  case  helium  and  iron  spectra 
[obtained  from  (eq.  4)]  are  multiplied  by: 

28.4  e-E/13.84  +  1.64  (12) 

for  E  _<  100  MeV/u.  The  resulting  spectra  are  employed  as  described  previously 
to  obtain  the  other  elemental  spectra,  i.e.  in  the  same  way  as  F^e  and  FFe 
were  used.  Again  this  worst-case  only  applies  for  E  <_  100  MeV/u,  above  that 
energy,  use  the  "pure"  galactic  cosmic-ray  spectra. 

In  addition  to  galactic  cosmic  rays,  some  particles  are  believed  to  be 
accelerated  in  the  interplanetary  medium.  The  most  important  of  these  is 
called  the  anomalous  component.  The  contribution  of  the  anomalous  component 
to  the  helium  spectrum  is  important  for  cosmic  ray  effects  on  micro¬ 
electronics.  We  recommend  that,  for  the  period  1582-1990,  the  cosmic  ray 
helium  spectrum  be  modified  as  follows: 

1.  Determine  the  maxima  values  of  the  cosmic-ray  spectra  from  Eq.  (4) 
using  the  He-max  and  He-min  constants  from  Table  1. 

2.  Modify  eq.  (4)  so  that  these  maximum  values  apply  for  all  energies 
below  the  energy  at  which  the  maxima  occurs,  ie.,  for  solar  minimum: 
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0.4  for 


E  <  195  MeV/u 


(13) 


f*He-m1n  *  j 

(  fHe-min  [from  eq.  (4)]  E  ■>  195  MeV/u 

Make  the  same  kind  of  modification,  f^-nuix  ^or  so1ar  maximum. 

3.  Combine  the  resulting  spectra  as  before  using  eqs.  (1-3). 

NOTE:  This  applies  only  to  He,  use  the  regular  He  spectra  eqs.  (1-6),  for 
obtaining  the  spectra  of  the  other  elements. 

Besides  helium,  the  anomalous  component  contributes  to  the  oxygen  and 
nitrogen  spectra  at  low  energies.  For  the  years  1982-1990,  these 
contributions  may  be  added  to  the  galactic  cosmic  ray  oxygen  and  nitrogen 
spectra  as  follows: 

For  oxygen,  use: 

f(E)  =  6  x  10“2exp[-( In  E-l. 79)2/0. 70]particles/m2ster. sec. MeV/u  (14) 

This  spectrum  crosses  over  the  galactic  spectrum  at  -  30  MeV.  The  two 
should  be  blended  at  that  point  with  eq.  (12)  replacing  the  galactic  spectrum 
at  lower  energies.  Similarly  for  nitrogen: 

f(E)  =  1.54  x  10_2exp[-(lnE-l. 79)2/0. 70]particles/m2ster. sec. MeV/u(15) 

Again,  this  blends  with  the  galactic  cosmic  ray  spectrum  at  about  30  MeV/u  and 
should  replace  it  below  this  energy. 

The  spectra  of  the  remaining  elements  are  unaffected  or  affected  at  too 
low  an  energy  to  matter. 

Solar  flare  particle  events  are  sporadic  occurrences  lasting  1-5  days. 
When  these  events  occur  they  can  be  the  dominant  cause  of  soft  errors.  For 
statistical  treatment,  they  are  broken  into  two  classes,  ordinary  (OR)  and 
anomalously  large  (AL).  The  probability  of  having  more  than  a  number  of 
events,  n,  in  a  time,  t,  is  given  by: 


P(n,t,N,T)  =  1  -  £  (i  +  N):(t/T)i/[i!N:(l  +  t/T)l+i+N]  (16) 

i=o 


where  T  and  t  are  in  years,  and  N  is  the  number  of  flares  that  have  occurred 
in  T  years. 

For  ordinary  events,  eq.  (16)  becomes: 

p0R  =  P(n,t,24,7)  for  (1981-1983) 

and  s  (17) 

PCR  =  P(n,t,6,8)  for  (1984-1988) 
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where  there  Is  a  probability  Por  of  having  more  than  n  ordinary  events  In  t 
years.  Similarly  for  anomalously  large  events: 

PAL  »  P(n,t,l,7)  (18) 

If  there  is  an  unacceptable  risk  of  an  AL  event  then  It  will  be  the  worst  case 
flare  for  the  mission. 

A  typical  OR  event-integrated  proton  differential  energy  spectrum  Is 
given  by, 

PtOR  *  3*3  x  ioS(e*£/20.2  +  307  e"E/3)  protons/m2ster.MeV  (19) 

where  E  Is  in  MeV. 

Ordinary  flares  come  In  a  broad  range  of  sizes.  With  a  ~  90  per  cent 
confidence  level,  an  ordinary  flare  spectrum  should  not  be  worse  than: 

fWOR  *  7.6  x  10^(e~8/30 .+  igs  e"8/^)  protons/m^ster.MeV  (20) 

where  E  Is  In  MeV. 


The  peak  proton  flux  differential  energy  spectrum  for  ordinary  events 
is,  typically: 

ftoR  =  1.95  x  io4(e-E/27.5  +  i73e-E/4)  protons/m2ster.sec.MeV  (21) 

where  E  is  in  MeY. 
and  no  worse  than: 

fW0R  a  1.71xl05(e-E/24.5+  63.6  e-E/4)  protons/m2ster.sec.MeV  (22) 

with  a  confidence  of  -  go  per  cent. 

Using  the  August  1972  flare  as  a  model  AL  event,  the  flare  integrated 
proton  differential  energy  spectrum  is: 

FAL  =  2.37  xlO11  exp[(30-E)/26.5]  protons/m2ster.MeV  (23) 

with  E  in  MeV  (King,  1974). 


The  peak  proton  flux  differential  energy  spectrum  is: 

I  9.3  x  10S  ( dP/dE)  exp(-P/0.10)  E  <  150  MeV 


fAL  a 


(24) 


\  1.76  x  105  (dP/dE)  P-9  E  >_  150  MeV 

in  protons/m2ster.sec.MeV,  where 


90 


P  -  C ( E/1000 )2  +  1.86  x  10-3e]1/2  (25) 

and  E  Is  In  MeY. 

The  composition  of  flare  particles  Is  also  highly  variable  from  flare 
to  flare.  Table  6  gives  the  composition  relative  to  hydrogen  for  the 
elements  through  nickel.  Both  mean  and  (90  per  cent  confidence  level)  worst 
cases  are  given.  To  obtain  the  spectrum  of  any  element  In  a  flare  just 
multiply  the  abundance  from  Table  6  by  the  appropriate  flare  proton  spectrum. 

TABLE  6.  Mean  and  Worst  Case  Compositions 


Mean  Case  Worst  Case 


H 

1 

1 

He 

2.2 

X 

IO"2 

3.3 

X 

H 

O 

9 

ro 

Li 

0 

0 

Be 

0 

0 

6 

0 

0 

C 

1.6 

X 

IO"4 

4.0 

X 

i 

o 

N 

3.8 

X 

IO"5 

1.1 

X 

1 

o 

rH 

0 

3.2 

X 

10~4 

1.0 

X 

10“3 

F 

0 

0 

He 

5.1 

X 

io-5 

l.S 

X 

ni- 

i 

o 

H 

Na 

1.6 

X 

LO 

1 

o 

1 — 1 

6.1 

X 

1C-6 

Mg 

4.8 

X 

l0-5 

1.8 

X 

1 

o 

r-( 

A1 

3.5 

X 

10~6 

1.4 

X 

10‘5 

Si  * 

3.8 

X 

10~5 

1.6 

X 

10-4 

Mean  Case 

Worst  Case 

p 

2.3 

x  10“7 

1.1 

x  IO"6 

s 

1.8 

x  10" 5 

8.4 

x  10"5 

Cl 

1.7 

x  10"7 

8  x 

10"7 

Ar 

3.9 

x  IO-6 

1.8 

x  10-5 

K 

1.3 

x  10"7 

6  x 

10"7 

Ca 

2.3 

x  10"6 

1  X 

10"5 

Sc 

0 

0 

Ti 

1  X 

io-7 

5  x 

10"7 

V 

0 

0 

Cr 

5.7 

x  10"7 

3.2 

x  10"6 

Mn 

4.2 

x  10"7 

2.3 

x  10"6 

Fe 

4.1 

x  10"5 

2.3 

*3- 

1 

o 

rH 

X 

Co 

1  X 

10-7 

5.5 

x  10"7 

Hi 

2.2 

x  10"6 

1.2 

x  10"5 

Z>28 

0 

0 

There  are  several  good  mathematical  models  for  the  trapped  proton 
environment.  We  recommend,  "AP-8  Trapped  Proton  Environment  for  Solar 
Maximum  and  Solar  Minimum"  by  Donald  M.  Sawyer  and  James  I.  Vette,  Report 
no.  NSSCS/WDC-A-R  and  S  7G-06,  Dec  1976,  NASA-Goddard,  Greenbelt,  Md. 


The  a  particle  (helium  nucleus)  flux  in  the  trapped  radiation  can  be 
scaled  from  the  proton  flux  by: 


a/p 


2.5x10-3 

8. 3xl0-3exp( -E/4. 15) 
0. 957exp( -E/1 .4) 


E  <  5  MeV/u 

5  <_  E  <  10  MeV/u  (25) 
E  >  10  MeV/u 
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We  suggest  that  this  formula  be  used  throughout  the  magnetosphere. 

Based  on  the  exceedingly  scanty  c;ata  available,  trapped  nuclei  heavier 
than  helium  are  probably  only  Important,  for  this  study,  In  the  outer  belt. 
For  L  >  3  (Mcllwaln’s  t  value),  and  suggest  that, 

Fcjho  *  3.8  x  10*  E"5.2r,ucie^/m2ster#sec,MeY/u  (27) 

with  E  In  MeV/u,  be  used  for  carbon  +  nitrogen  +  oxygen. 

The  modulation  of  cosmic  ray  spectra  by  the  earth's  magnetic  field 
requires  a  more  thorough  treatment  than  can  be  offered  here,  but  some 
guidance  will  be  provided.  The  geomagnetic  cutoff  Is  a  value  of  magnetic 
rigidity  below  which  cosmic  rays  will  not  reach  a  specified  point  in  the 
magnetosphere  from  a  specified  direction.  The  magnetic  rigidity,  P,  in 
Gev/ec  may  be  computed  from  the  particles'  energy  using: 

P  =  (E/1000)2  +  1.86  x  10“3e]1/2  (28) 

where  E  is  in  MeV/u  and  A  and  2  are  the  atomic  mass  and  charge  of  the 
nucleus  in  question. 

The  cutoff  is  most  simply  computed  with: 

Pc  =  15 . 96/L2 -  005  (29) 

where  L  is  Me II wain's  L  parameter 

Detailed  calculations  of  the  cutoff  are  available  from  Shea  and  Smart 
( 1975 ) .  Transmittance  functions  for  satellite  orbits  may  be  computed  using 
the  techniques  described  in  Heinrich  and  Spill,  (1S79). 

The  transmittance  functions  of  Heinrich  and  Spill  are  useful  in 
moGulating  the  cosmic  ray  spectra.  Some  thought  must  be  given  to  their  use 
on  solar  flare  spectra  because  the  flare  particle  intensity  changes  on  a 
time  scale  comparable  to  or  shorter  than  an  orbital  period.  Also  the 
geomagnetic  cutoff  is  suppressed  to  some  extent  during  a  flare.  We 
recommend  that  the  geomagnetic  cutoff  during  a  flare,  Pr,  be  computed  from 
the  "quiet  time"  cutoff  P0  using: 

aP/P0  =  0.54  exp(-P0/2.S)  (30) 

and 

PF  =  po  -  *p  (31) 

where  Pp,  PQ  and  aP  are  in  GeV/ec. 
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